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Abstract. Let R be a real closed field and D C R an ordered do- 
main. We give an algorithm that takes as input a polynomial Q C 
D[Ai, . . . , Xk], and computes a description of a roadmap of the set of ze- 
ros, Zer(Q,R fe ), of Q in R fc . The complexity of the algorithm, measured 
by the number of arithmetic operations in the domain D, is bounded 
by d ol - k ^\ where d = deg(Q) > 2. As a consequence, there exist 
algorithms for computing the number of semi-algebraically connected 
components of a real algebraic set, Zer(Q, R fe ), whose complexity is also 
bounded by d ol - k ^\ where d — deg(Q) > 2. The best previously known 
algorithm for constructing a roadmap of a real algebraic subset of R fe 
defined by a polynomial of degree d had complexity d°' fc '. 



1. Introduction 

The problem of designing efficient algorithms for deciding whether two 
points belong to the same semi-algebraically connected component of a semi- 
algebraic set, as well as counting the number of semi-algebraically connected 
components of a given semi-algebraic set ScR fc where R is a real closed 
field (for example the field of real numbers), is a very important problem in 
algorithmic semi-algebraic geometry. 

The first algorithm for solving this problem [12] was based on the tech- 
nique of cylindrical algebraic decomposition [TJ [2], and consequently had 
doubly exponential complexity. 

Algorithms with singly exponential complexity were given later in a series 
of papers EJ QUI QU CQ . 

They are all based on a geometric idea introduced by Canny, the construc- 
tion of an one-dimensional semi-algebraic subset of the given semi-algebraic 
set S, called a roadmap of S, which has the property that it is non-empty 
and semi-algebraically connected inside every semi-algebraically connected 
component of S. 
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In the papers mentioned above, the construction of a roadmap of a semi- 
algebraic set S depends on recursive calls to itself on several (in fact, singly 
exponentially many) (k — 1) dimensional slices of S, each obtained by fixing 
the first coordinate. For constructing the roadmap of a real algebraic variety 
defined by a polynomial Q C R[-X~i, . . . , X^] with deg(Q) < d, this technique 
gave an algorithm with complexity d°^ k ) . The exponent in the complexity, 
0(/c 2 ), is due to the fact that the depth of the recursion in these algorithms 
could be as large as k. This exponent is not satisfactory since the total 
number of semi-algebraically connected components is (0(d)) k and so there 
is room for trying to improve it. However, this has turned out to be a rather 
difficult problem with no progress till very recently. 

A new construction for computing roadmaps, with an improved recursive 
scheme of baby step - giant step type, has been proposed, and applied suc- 
cessfully in the case of smooth real algebraic hypersurfaces in [9]. In this 
new recursive scheme, the dimension drops by \fk in each recursive call. As 
a result, the depth of the recursive calls in this new algorithm is at most \fk, 
and consequently the algorithm has a complexity of rf°( fcv ^). The proof of 
correctness of the algorithm in [9] depends on certain results from commuta- 
tive algebra and complex algebraic geometry, in order to prove smoothness 
of polar varieties corresponding to generic projections of a non-singular hy- 
persurface. Choosing generic coordinates in the algorithm is necessary since 
the non-singularity of polar varieties does not hold for all projections, but 
only for a Zariski-dense set of projections. This is an important restriction, 
since there is no known method for making such a choice of generic coordi- 
nates deterministically within this improved complexity bound. As a result, 
the authors obtain a randomized (rather than a deterministic) algorithm for 
computing roadmaps: there might be cases where the algorithm terminates 
and gives a wrong result. 

In contrast to these techniques which depend on complex algebraic ge- 
ometry, the algorithm for constructing roadmaps described in [2] depend 
mostly on arguments which are semi-algebraic in nature. The greater flexi- 
bility of semi-algebraic geometry (as opposed to complex geometry) makes 
it possible to avoid genericity requirements for coordinates. More precisely, 
we apply the technique used in [2] to make an infinitesimal deformation 
of the given variety so that the original coordinates are good. Since the 
infinitesimal deformation uses only one infinitesimal, it does not affect the 
asymptotic complexity class of the algorithm. 

The goal of this paper is to obtain a deterministic algorithm for computing 
the roadmap of a general algebraic set, combining a baby step - giant step 
recursive scheme similar to that used in [9] and extending techniques coming 
from [2]. 



We start by recalling the precise definition of what is meant by a roadmap. 
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Definition 1.1. Let S C R fc be a semi-algebraic set. A roadmap for S is 
a semi-algebraic set RM(S') of dimension at most one contained in S which 
satisfies the following roadmap conditions: 

(1) RMi For every semi-algebraically connected component C of S, C n 
RM(jS') is semi-algebraically connected. 

(2) RM2 For every x € R and for every semi-algebraically connected 
component D of S x , DnRM(S') 7^ 0, where we denote by S x the set 
S n ir^ 1 (x) for x G R, and 7Ti : R fc — > R the projection map onto the 
first coordinate. 

Let M C R fc be a finite set of points. A roadmap for (S,A4) is a semi- 
algebraic set RM(«S, M) such that RM(6 r , Ai) is a roadmap of S and M C 
RM(S,M). 

The main result of the paper is the following theorem. The notion of real 
univariate representations used in the following statements is explained in 
Section |H 

Theorem 1.2. Let Z C R fc be an algebraic set defined as the set of zeros of 
a polynomial of degree at most d > 2 in k variables with coefficients in an 
ordered domain D contained in a real closed field R. 

a) There exists an algorithm for constructing a roadmap for Z using 
^O(fevfe) arithmetic operations in D. 

b) Moreover, there exists an algorithm that given a finite set of points 
Mo C Z , with cardinality 5, and described by real univariate repre- 
sentations of degree at most d°^ k \ constructs a roadmap for (Z, A4o) 
using 5°^ jO(k^/k) arithmetic operations in D. 

The following corollary is an immediate consequence of b). 

Corollary 1.3. Let Z C R fc be an algebraic set defined as the set of zeros 
of a polynomial of degree at most d > 2 in k variables with coefficients in 
an ordered domain D contained in a real closed field R. 

a) There exists an algorithm for counting the number of semi-algebraically 
connected components of Z which uses d°^ k ^ arithmetic operations 
in D. 

b) There exists an algorithm for deciding whether two given points, de- 
scribed by real univariate representations of degree at most dP^ k ' be- 
long to the same semi-algebraically connected component of Z which 
uses d°( fcv/ k) arithmetic operations in D. 

Remark 1.4. We can always suppose without loss of generality that the zero 
set of a family of polynomials of degree at most d is defined by one single 
polynomial of degree at most 2d by replacing the input polynomials by their 
sum of squares. 

Remark 1.5. Even if the input is a polynomial with coefficients in the field of 
real numbers, the deformation techniques by infinitesimal elements we use 
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make it necessary to perform computations on polynomials with coefficients 
in some non-archimedean real closed field. This is the reason why general 
real closed fields provide a natural framework for our work. 

2. Outline 

We outline below the classical construction of a roadmap RM(Zer(£j), R fc )) 
for a bounded algebraic set Zev(Q,K k ), defined as the zero set of a polyno- 
mial Q inside R fc . The geometric ideas yielding this construction are due 
to Canny. The description below is similar to the one in [21 Chapter 15, 
Section 15.2]. 

A key ingredient of the algorithm is the construction of a particular finite 
set of points intersecting every semi-algebraically connected component of 
Zer(Q,R fc ). In the case of a bounded and non-singular real algebraic set 
in R fe (in the generic case), these points are nothing but the set of critical 
points of the projection to the Xi-coordinate on Zer(Q, R fc ). In more general 
situations, the points we consider are called X\ -pseudo-critical points, since 
they are obtained as limits of the critical points of the projection to the 
Xi-coordinate of a bounded nonsingular algebraic hypersurface defined by a 
particular infinitesimal deformation of the polynomial Q. Their projections 
on the Xi-axis are called pseudo-critical values. 

We first construct the "silhouette" which is the set of -XVpseudo-critical 
points on Zer(Q, R fc ) along the Xi-axis by following continuously, as x varies 
on the Xi-axis, the X2-pseudo-critical points on Zer(Q, R k ) x . This results 
in curves and their endpoints on Zer(<5,R fc ). The curves are continuous 
semi-algebraic curves parametrized by open intervals on the Xi-axis and 
their endpoints are points of Zer(Q,R fc ) above the corresponding endpoints 
of the open intervals. Since these curves and their endpoints include for 
every x £ R the X2-pseudo-critical points of Zer(Q, R fc ) x , they meet every 
semi-algebraically connected component of Zer(Q,R k ) x . Thus, the set of 
curves and their endpoints, already satisfy RM2. However, it is clear that 
this set might not be semi-algebraically connected in a semi-algebraically 
connected component and so RMi might not be satisfied. 

In order to ensure property RMi we need to add more curves to the 
roadmap. For this purpose, we define the set of distinguished values T> as 
the union of the X\ -pseudo-critical values, and the first coordinates of the 
endpoints of the curves described in the previous paragraph. A distinguished 
hyperplane is an hyperplane defined by X\ = v, where v is a distinguished 
value. The input points, the endpoints of the curves, and the intersections of 
the curves with the distinguished hyperplanes define the set of distinguished 
points, M.. 

Let the distinguished values be V\ < . . . < ujy. Note that amongst these 
are the Xi-pseudo-critical values. Above each interval (vi, fj+i) we have con- 
structed a collection of curves Ci meeting every semi-algebraically connected 
component of Zer(Q,R fc )„ for every v € (v{, Above each distinguished 
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value Vi we have a set of distinguished points Mi- Each curve in Ci has an 
endpoint in Mi and another in Ni+\. Moreover, the union of the A/i contains 
A/". We denote by C the union of the Ci . 

The following key connectivity result is proved in [31 Lemma 15.9]. 

Proposition 2.1. Let TZ = C U Zer(Q,R k ) v . If P is a semi- algebraically 
connected component of Zer(Q, R fc ), then TZ n P is semi- algebraically con- 
nected. 

Thus, in order to construct a roadmap of Zer(Q,R fc ) it suffices to re- 
peat the same construction in each distinguished hyperplane Hi defined by 
X\ = Vi with input Q(vi,X2, ■ ■ ■ and the distinguished points in M. Vi 
by making recursive calls to the algorithm. The following proposition is 
proved in [21 Proposition 15.7]. 

Proposition 2.2. The semi- algebraic set RM(Zer(Q, R fc ), M) obtained by 
this construction is a roadmap for Zer(Q,R fc ) containing M. 

To summarize, classical roadmap algorithms based on Canny's construc- 
tion proceed by first considering the "silhouette" , consisting of curves in the 
Xi-direction, and then making recursive calls to the same algorithm at cer- 
tain hyperplane sections of Zer(Q, R fc ), so that the dimension of the ambient 
space drops by 1 at each recursive call. 

The main difference between classical roadmap algorithms and the algo- 
rithms described in [9] and in the current paper is that instead of considering 
curves in the Xi-direction and making recursive calls to the same algorithm 
at certain hyperplane sections of Zer(Q,R fc ) corresponding to special values 
of X±, so that the dimension of the ambient space drops by 1, we consider 
a p-dimensional subset W of Zer(Q,R fc ) where 1 < p < k, and make re- 
cursive calls at certain (k — p)-dimensional fibers of Zer(Q,R fc ), so that the 
dimension of the ambient space drops by p. 

The main topological result, generalizing Proposition ^. 1\ is that the semi- 
algebraic set which is the union of W and these fibers is semi-algebraically 
connected. This is proved in Section El in a special case. Thus, in order 
to produce a roadmap of Zer(Q,R fe ) it suffices to compute a roadmap of 
W passing through an appropriate set of points, and the roadmap of the 
corresponding fibers in a (k — p)-dimensional ambient space, using recursive 
calls. 

The fact that in the new algorithm we are fixing a whole block of p 
variables at a time necessitates introducing a new kind of algebraic repre- 
sentation which we call "real block representation". This notion is defined 
in Section HI where we also explain how to represent curves. 

In Section (U the roadmap of W is computed by an algorithm directly 
adapted from [2j Algorithm 15.3] which makes use of the fact that W is 
low dimensional, in a special case. The general case, requiring the use of a 
deformation technique and a limit process, is described in Section El 
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Finally, we obtain in Section [7] a baby step - giant step roadmap algorithm 
for a general algebraic set. We prove its correctness, as well as the improved 
complexity bound. 

The algorithm for computing efficiently limits of curve segments is quite 
technical. Since, this technicality can obscure the ideas behind the main 
algorithm, for the sake of readability we have postponed the details behind 
taking limits of curves to a separate section (Section [8|). 

Throughout the paper, we use as a basic reference [2]. We cite [3] instead 
when the precise statements needed in the paper are not included in [2] . 

3. Connectivity results 

In this section we prove a topological result about connectivity which will 
be used in proving the correctness of our algorithm later. The statement 
of the result, as well as the main ideas of the proof, is influenced by [9J 
Theorem 14]. It is a direct generalization of Proposition 12.11 to the case of 
projection onto more than one variable. 

We denote by R a real closed field. 

Notation 3.1. For 1 < q < p < k, we denote by 7rr 5p i : R fc ->■ W~ g+1 the 
projection 

\X\ , . . . , Xk ) I }■ (yXq , . . . , Xp) . 

In case p = q we will denote by tt p the projection 7rrppi. For 1 < q < p < k, 
we denote by flVp] : R fc = R fc — > W~ q the projection 

(Xl, . . . , X/j) I y (Xq+i, . . . , Xp). 

For any semi-algebraic subset S C K k , and T C R p , we denote by St the 
semi-algebraic set ttZXT) n S, and S y rather than S{ y y, for y G R p . We 
also denote S <a and S< a rather than 5(_ 00)0 ) and S , (_ 00)(1 ], for a G R. 

We denote as before by Zer(Q, R fc ) the algebraic set of zeros of a poly- 
nomial Q G R[Xi, . . . ,Xf-] inside R fc . Note that this does not imply that 
the dimension of Zer(<5,R fc ) is A; — 1. In fact, over any real closed field, 
algebraic sets defined by one equation coincide with general algebraic sets 
since replacing several equations by their sum of squares does not modify 
the zero set. A Q-singular point is a point x such that 

«,)-^ W -...-^(x)-a 

Note that this is an algebraic property related to the equation Q rather than 
a geometric property of the underlying set Zev(Q,K k ): two equations can 
define the same algebraic set but have a different set of singular points. 
Similarly a Q-critical point of 7Ti is a point x such that 
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To simplify notations, when there will be no ambiguity on Q, we will simply 
refer to singular /critical points. 

In this paper, we will be using constantly the notion of semi-algebraically 
connected components of a semi-algebraic set [2j Section 5.2]. Note that, in 
particular, a semi-algebraically connected component is always non-empty 
by definition [3j Theorem 5.21]. 

Property 3.2. We now consider a tuple 

(V,p,W,(Mi)i<i<2,(Vi)i<i<2) 
with the following properties 

(1) V C R k is the union of certain bounded semi-algebraically connected 
components of an algebraic set Zer(Q,R fc ) C R fc , such that the Q- 
singular points of V, as well as the Q-critical points of the map tti 
on V form the finite set M.\ C V, and T>\ = iri(Mi); 

(2) W C V is a closed semi-algebraic set of dimension p, 1 < p < k, such 
that for each y G W, W y is a finite set of points having non-empty 
intersection with every semi-algebraically connected component of 

v y ; 

(3) M.2 C V is a finite subset such that the intersection of At 2 with every 
semi-algebraically connected component of W a is non-empty, for a £ 
T^2 — 7r i(-A / ^2)- Moreover for every interval [a,b] and c G [0,6] with 
{c} D 2?2n [a, b], if D is a semi-algebraically connected component of 
VF[ a w, then D c is a semi-algebraically connected component of W c . 

A tuple 

(V,p,W,{Mi)l<i< 2 ,(T> i ) 1 < i <2) 

satisfies Property 13.21 if it satisfies the above properties (1) to (3). 

We state now the main result of this section. It generalizes Proposition 
12.11 as well as [9[ Theorem 14], in the special case of Property 13.21 

Proposition 3.3. Let 

(V,p,W,(M)l<i<2,(A)l<i<2) 

satisfy Property \3. 21 

M = T[ij,](-Mi),M = 7r[i iP ](A^2), 

and 

S = WUVviuATa- 

For every semi-algebraically connected component C of V , C f] S is non- 
empty and semi-algebraically connected. 

Remark 3.4. In order to understand the situation, the following example of 
a tuple satisfying Property 13.21 can be useful 
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(1) the torus V CR 3 denned as the set of zeros of the equation 

q = nxf + ( 12X \ + 3 5X3 ) 2 - (*? + x " + *3 + s) 2 , 

([1], page 40, figure 2.5), the four critical points M.\ C V, of the 
map 7Ti restricted to V and T>\ = 7i"i(A4i); 

(2) the silhouette W C V defined by 

(3) the six critical values T>2 C R of the map n\ restricted to the sil- 
houette W, and the intersection M2 of the corresponding six fibers 
with the silhouette W. 

The tuple 



(y,l,W,(Mi)i<i<2,(Vi] 



Ki<2, 



satisfies Property | 

Finally, S is the union of the silhouette and the intersection of the torus 
with the six curves which are the fibers of V at the distinguished values T>2- 



The end of this section is devoted to prove Proposition 13.31 We need pre- 
liminaries about non-archimedean extensions of the base real closed field R. 

Remark 3.5. A typical non-archimedean extension of R is the field R(e) 
of algebraic Puiseux series with coefficients in R, which coincide with the 
germs of semi- algebraic continuous functions (see [21 Chapter 2, Section 6 
and Chapter 3, Section 3]). An element x G R(e) is bounded over R if \x\ < r 
for some < r G R. The subring R(e}& of elements of R(e) bounded over 
R consists of the Puiseux series with non-negative exponents. We denote 
by lim e the ring homomorphism from R(e}& to R which maps X^eN *^ 9 
to ao- So, the mapping lim £ simply replaces e by in a bounded Puiseux 
series. Given S C R(e) fc , we denote by lim e (5) C R fc the image by lim e of 
the elements of S whose coordinates are bounded over R. 

More generally, let R' be a real closed field extension of R. If S C R fc is a 
semi-algebraic set, defined by a boolean formula $ with coefficients in R, we 
denote by Ext (5, R') the extension of S to R', i.e. the semi-algebraic subset 
of R /fc defined by <3?. The first property of Ext(S l , R') is that it is well defined, 
i.e. independent on the formula $ describing S [21 Proposition 2.87]. Many 
properties of S can be transferred to Ext (5, R'): for example S is non-empty 
if and only if Ext (5, R') is non-empty, S is semi-algebraically connected if 
and only if Ext (5, R') is semi-algebraically connected [21 Proposition 5.24]. 

Moreover, if Property 13.21 (2) holds for V, W, i.e. for every y € R p , 
W y is a finite set of points having non-empty intersection with every semi- 
algebraically connected component of V y , then Property 13.21 (2) holds for 
Ext(V,R'),Ext(W,R'), i.e for each y' G R /p , Ext(W,R'V is a finite set 
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of points having non-empty intersection with every semi-algebraically con- 
nected component of Ext(V, R')?/- Indeed, by Hardt's semi-algebraic trivial- 
ity theorem [2J Theorem 5.45], one can find a finite partition of W in semi- 
algebraic sets Tj, i = 1, . . . , r, a finite partition of Vjv into semi-algebraic 
sets Sij and an integer nj > such that Sij is semi-algebraically home- 
omorphic to T x (Sij) yi for some y, G Tj, and for all y G Tj, the semi- 
algebraically connected components of V y are (Sij) y and Wy has raj points. 
By Tarski-Seidenberg's transfer principle [21 Theorem 2.80], Ext(Sjj, R') is 
semi-algebraically homeomorphic to Ext(T,R') x Ext(<Sjj, R').^ , and for all 
y' G R /p , there exists i such that y' G Ext(T, R'), the sets Ext(S'jj, R')i/' are 
the semi-algebraically connected components of Ext(V,R')j / / and the inter- 
section of Ext(W, R')y/ and Ext(S , j J , R') y / has exactly rij points. 

We now prove a few preliminary results about V defined as the union 
of certain bounded semi-algebraically connected components of an algebraic 
set Zer(Q, R fc ) C R fc , supposing that the set Ai± of points which are singular 
points or critical points of 717 on Zer(Q, R fc ) inside V is finite. 

In this paper a semi- algebraic path is a semi-algebraic continuous func- 
tion 7 from a closed interval [a, b] C R to R*. Note that a semi-algebraic 
set is semi-algebraically connected if and only it is semi-algebraically path 
connected [21 Theorem 5.23]. 

Lemma 3.6. Suppose that b G" T>\ = Let C be a semi-algebraically 

connected component ofV<b- If a < b and (a, b] nV\ is empty, then C< a is 
semi-algebraically connected. 

Proof. Let x and y be two points of C< a and 7 : [0, 1] — > C be a semi- 
algebraic path connecting x to y inside C. We want to prove that there is a 
semi-algebraic path connecting x to y inside C< a . 

If Im(7) C C< a there is nothing to prove. 

If Im( 7 ) <£ C~< a , 

3c G R Va < d < c Tm(i) n Zer(Q) d / 0. 
Let e be a positive infinitesimal. Then 

Ext( 7 ([0, 1]), R(e» n Zer(Q, R(e) k ) a+£ / 
using [21 Proposition 3.17] . Since 

{u G [0,1] C R(e) I Ext(7,R(e))(«) G Zer(Q,R( £ ) fc ) <0+£ } 

and 

{n G [0,1] C R(e) I Ext( 7 ,R(e))(n) G Zer(Q, R(e) k ) [a+£M } 
are semi-algebraic subsets of [0, 1] C R(e) there exists by [21 Corollary 
2.79] a finite partition *p of [0, 1] C R(e) such that for each open inter- 
val (u,v) of qj, Ext(7,R(e))(« , f ) is either contained in Zer(Q, )<a+e; 
or in Zer(Q,R(e) fc )[ a+£i6 ], with -y(u) and 7(u) in C a+£ . 

If Ext (7, R(e))(«, u) is contained in Zer(Q, R(e) fc )[ a+e> b], we can replace 7 
by a semi-algebraic path 7? & , connecting 7(14) to 7(f) inside C a + e . Note that 
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there is no Q-critical point of m in Ext(V, R(e))r 0+e w and Ext(V, R(£})r a +E,&] 
contains no Q-singular point by [21 Proposition 3.17] while Ext(V, R(e)) C 
Zer(Q,R(e) fc ) by Proposition 2.87] . 

By [31 Proposition 15.1 b] if D is a semi-algebraically connected compo- 
nent of Ext(V, R(e))[ a+e u , D a+£ is a semi-algebraically connected compo- 
nent of Ext(V,R(e)) a+£ .' 

Construct a semi-algebraic path 7' from x to x' inside C< a + e , obtained by 
concatenating pieces of 7 inside Zer(Q, R(e) fc ) <a+£ and the paths *ft uv \ con- 
necting 7(u) to 7(f) for (u, v) such that Ext (7, R(e))(u, v) C Ext(V, R(e))[ a+e 6 ] 
Note that such a semi-algebraically connected path 7' is closed and bounded. 
Applying [5J Proposition 12.43], lim e (7'([0, 1])) is semi-algebraically con- 
nected, contains x and x' and is contained in lim e (C< a+e ) = C< a . This is 
enough to prove the lemma. □ 

We continue to suppose that M\ is finite. 

Lemma 3.7. Let C be a semi-algebraically connected component of V<b ■ 
such that C CiVb is not empty. 

(1) If dim(C) = 0, C is a point contained in M\. 

(2) If dim(C) 7^ 0, C<fe is non-empty. Let B\, . . . ,B r be the semi- 
algebraically connected components of C<&. Then, 

(a) for each i, 1 < i < r, L\nMi / 0; 

(b) if there exist < i < j < r such that Bi n Bj ^ 0, then 
B~nB~ C Mi; 

(c) yJ' i=1 Bi = C , and hence U r i=l Bi is semi-algebraically connected. 

Proof. Part 1 follows immediately from [31 Proposition 7.3]. Let us prove 
Part 2: since M\ is finite, there is a non-singular point x G C which is non- 
critical for 7Ti on V . Let T X V denote the tangent space to 7 at 1. So T X V 
is not orthogonal to the X\ axis, and the semi-algebraic implicit function 
theorem [21 Theorem 3.25] implies that C<b is non-empty. 

Part 2) a) and 2 b) are immediate consequences of Proposition 7.3 in [3]. 

We prove 2) c). Clearly, U r i=1 B~ C C. Suppose that x E C \ U r i=1 ~B~. 
For r > and small enough, Bk(x,r) n C<fc = (where Bk(x,r) is the k- 
dimensional open ball of center x and radius r). Note that tti(x) = b, since 
otherwise x belongs to C<b, and thus to one of the B^s. 

Applying [31 Proposition 7.3], we deduce from the fact that Bk(x,r) n 
C<f, = Bk(x, r)<6 n C = that x is either a Q-singular point, or a Q-critical 
point of 7Ti on V. In other words x G A4i. But since by assumption .Mi is 
finite, this implies that C\J r i=l Bi is a finite set. Since C is semi-algebraically 
connected and of positive dimension, C \ yj r i=1 Bi must be empty. □ 

Notation 3.8. If S C R fc is semi-algebraic set and x G S, then we denote 
by C(S,x) the semi-algebraically connected component of S containing x. 

We are now ready to prove Proposition 13.31 
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Proof of Proposition \3.3[ For a in R, we say that property P(a) holds if: for 
any semi-algebraically connected component C of V< a , C n S is non-empty 
and semi-algebraically connected. 

We prove that for all a in R, P(a) holds; taking o > max^gy 7Ti(x) suffices 
to prove the proposition since V is bounded. 

Let V = Vi U V 2 = m(Mi U M 2 ) (see Property E2} • 

The proof uses two intermediate results: 
Step 1: For every aeP, P(o) implies P(6) for all b G R with (a, 6] nD = 0. 
Step 2: For every 6 G 2?, if P(a) holds for all a < b, then P(fe) holds. 

Since for a < min^y 7Ti(x), property P(a) holds vacuously, and the com- 
bination of these two results gives by an easy induction P(a) for all a in R, 
thereby proving the proposition. 

We now prove the two steps. 
Step 1. We suppose that a G T>, P(a) holds, take b G R, a < b with 
(a,b]nV = $ and prove that P(b) holds. Let C be a semi-algebraically con- 
nected component of V<b. We have to prove that CnS is semi-algebraically 
connected. 

Since (a, b] D X> = 0, it follows that (-Mi)(a,b] = ®i an d C< a is a semi- 
algebraically connected component of V< a using Lemma 13.61 So, using 
property P(a), we see that C< a n 5 is non-empty and semi-algebraically 
connected. 

If C< a n 5 = C n 5, there is nothing to prove. Otherwise, let x G C n 5 
such that X C< a . We prove that x can be semi-algebraically connected to 
a point in C< a n5 by a semi-algebraic path in C n 5, which is enough to 
prove that C n 5 is semi-algebraically connected. 

Since 7i"i(x) G (a, 6] and (a, 6] n 2? = 0, tti(x) P and x KviuA/i- So, 
from x G S, we get x G W. We note that C(Wt a ^, x) C C. By Property 13.21 
(3) applied to C(W[ a ^j, x) we have that a G 7Ti(C(W[ a ^, x)) and C(V^[ 0) &],:c) a 
is non-empty. Hence there exists a semi-algebraic path connecting x to a 
point in C(W[ a ^,x) a inside C{W\ a u,x). Since C{Wi a ^ux) C W C 5 and 
C(Wr a M, x) C C, if follows that C(Wy a ^x) C Cn5 and we are done. 
Step 2. We suppose that 6 G V, and P(a) holds for all a < b, and prove 
that P(6) holds. 

Let C be a semi-algebraically connected component of y<&. If Cf, = 
there is nothing to prove. Suppose that Cb is non-empty; we have to prove 
that C n S is semi-algebraically connected. 

If dim(C) = 0, C is a point, belonging to .Mi C 5 by Lemma [3.71 So 
C n 5 is semi-algebraically connected. 

Hence, we can assume that dim(C) > 0, so that C<& is non-empty by 
Lemma 13.71 

Our aim is to prove that C D S is semi-algebraically connected. We do 
this in two steps. We prove the following statements: 

(a) If B is a semi-algebraically connected component of C<6, then BnS 
is non-empty and semi-algebraically connected, and 



12 



S. BASU, M-F. ROY, M. SAFEY EL DIN, AND E. SCHOST 



(b) and, using (a) C nS is semi-algebraically connected. 

Proof of (a) We prove that if B is a semi-algebraically connected compo- 
nent of V<b, then B n S is non-empty and semi-algebraically connected. 

Since B contains a point of Af i it follows that B n S is not empty. 

Note that HBDS = BnS, then there exists a with 

max({7Ti(x) | x G B n 5}) < a < b, 

with £> n5 = (5n5)< a and B< a semi-algebraically connected using Lemma 
13.61 So B n S is semi-algebraically connected since P(a) holds. 

We now suppose that (B \ B) n 5 is non-empty. Taking x G (B\B)nS, 
we are going to show that x can be connected to a point z in n S by a 
semi-algebraic path 7 inside Bn5. Notice that 7Ti(x) = b. 

We first prove that we can assume without loss of generality that x G W. 
Otherwise, since x G S and S = W U V/ViUA/^ we must have that x G V y 
with y = vr[ lp ](x), and V y C 5. Let A = C(V y n B,x). We now prove 
that A n Wy 7^ 0. Using the curve section lemma choose a semi-algebraic 
path 7 : [0,e] — > Ext(J3,R(e}) such that 7(0) = x, lim e 7(e) = x and 
7((0,e]) C Ext(S,R(e)). Let y e = vr^^e)) and 

A B = C(Ext(J3,R(e)) v „7(e)). 

Note that x G lim e A £ C A. 

By Remark 13.51 Ext(-B,R(e)) is a semi-algebraically connected compo- 
nent of Ext(V< a , R(e)) which implies that A £ is a semi-algebraically con- 
nected component of Ext(V, R(e)) 2/e . By Property 13.21 (2) and Remark 13.5} 
Ext(W / ,R(e)) ?/E n A £ ^ 0. Then, since Ext(W,R(e)) ye n A £ is bounded over 
R, lim e (Ext(W, R{e)) y£ fl A e ) is a non-empty subset of W y n A 

Now connect x to a point in x' G W y by a semi-algebraic path whose 
image is contained in A C B y C (B \ B) P\ S. Thus, replacing x by x' if 
necessary we can assume that x £ W as announced. 

There are four cases, namely 

(1) x G Mi UM2; 

(2) x^MiUM 2 and C(W b , x)Ji B; 

(3) i^MiU Al 2 , C(W 6 , x) C 5 and 6 G £> 2 ; 

(4) x All U M 2 , C(W b , x) C B and 6 £> 2 ; 

that we consider now. 

(1) x G All UX 2 : 

Define y = 7rri jP i(x) G R p , and note that V y C S. Since x G B, 

and 5 is bounded, y G 7rn p i(-B) = 7Tm iP ](B). Now let e > be 
an infinitesimal. By applying the curve selection lemma to the 
set B and x G B, and then projecting to W using n\i jP ] we ob- 
tain that there exists y e G R(^) p infinitesimally close to y with 
TTi(jfe) < 7Ti(y), andx G lim e Ext(V, R(e)) 2/s . Let x £ G Ext(V, R(e)) ye 
be such that 7rhfc](limx E ) = x. Moreover, by Property 13.21 (2) 
and Remark 13.51 we have that Ext(W, H{e)) ys is non-empty and 
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meets every semi-algebraically component of Ext(V, R(e)) 3/e . . Let 
x' E G Ext(W,R(e))j /£ n C(Ext(fi, R(e)) ys , x 6 ), and x' = 7T[ 1)fc ] (lim E x E ). 
Since ttm w (lim x £ ) = x and lim E C(Ext(£>, K(e)) ys , x E ) is semi-algebraic- 
ally connected, 

limC(Ext( J B,R(e)) 2/e ,x £ ) cC(B y ,x). 

Now choose a semi-algebraic path 71 connecting x to x' inside C(B y , x) 
(and hence inside S since C(B y ,x) C V y C. S), and a semi-algebraic 
path 72(e) joining x' to x £ inside Ext(W, R(e)). The concatenation 
of 71 , 72 (e) gives a semi-algebraic path 7 having the required prop- 
erty, after replacing e in 72(e) by a small enough positive element of 
R. 

(2) x Mi U M 2 and C(W b , x) <£ B: _ 

There exists x' G C(Wft,a;), x' ^ B and a semi-algebraic path 7 : 
[0, 1] -»• C(W b ,x), with 7(0) = 1,7(1) = x'. Since x' B, it follows 
from Lemma ETTJ 2) that for i x = max{0 < £ < 1 | 7(4) G B}, 
7(^1) G A^i. We can now connect a/ to a point in B n 5 by a 
semi- algebraic path inside B (~)S using (1). 

(3) x # Mi U M 2 , C(W b , x) C 5 and b G £> 2 : 

Since 6 G £> 2 by Property El (2) there exists x' G C(TU 6 ,x) n A4 2 . 
Thus, there exists a semi-algebraic path connecting x to x' G M 2 
with image contained in B (~)W <Z B (~) S . We can now connect x' to 
a point in B n 5 by a semi-algebraic path inside 6n5 using (1). 

(4) xgMiU M 2 , C(W b , x) C B and 6 £> 2 : 

Since 6 £>2> for all a < b such that [a, 6] n V 2 = 0, C(Wr a M,x)& = 
C(W b ,x) &ndC{W [aM ,x) a / by Property (3). Let x' G C(W [aM ,x) a . 
We can choose a semi-algebraic path 7 : [0,1] — > C(W\ a w,x) with 
7(0) = 1,7(1) = x'. Let ti = max{0 < t < 1 | 7('i) G 
Then, either 7(ti) G Mi and we can connect 7(^1) to a point in 
B n 5 by a semi-algebraic path inside B OS using (1). Otherwise, 
by Lemma 13771 (2 b), for all small enough r > 0, 6^(7(41), r) n C<b is 
non-empty and contained in B. Then, there exists t 2 G (t±, 1] such 
that 7(^2) £B(lWcBnS, and the semi-algebraic path tI [o,* 2 ] 
gives us the required path in this case. 

Taking x and x 1 in B S, they can be connected to points z and z' in 
B n 5 by semi- algebraic path 7 and 7' inside Bfl5 such that, without loss 
of generality, tti(z) = 717 (z')) = a. Using P(a), we conclude that P(6) holds. 
Proof of (b) We have to prove that C n S is semi-algebraically connected. 

Let x and be in C n 5. We prove that it is possible to connect them 
by a semi-algebraic path inside C n <S. 

Since we suppose that dim(C) > 0, C <b is non-empty by Lemma 13.71 (2). 
Using Lemma [3771 (2.c). let Bi (resp. Bj) be a semi-algebraically connected 
component of C <b such that x G Bi (resp. x' G J3j). 
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If i = j, x and x' both lie in B^nS which is semi-algebraically connected 
by (a). Hence, they can be connected by a semi-algebraically connected 
path mBiDScCnS. 

So let us suppose that i ^ j. Note that: 

• by Lemma 13771 (2. a), Bi n M.\ and Bj n M.\ are not empty, 

• by (a) BidS and Bj n S are semi-algebraically connected, 

• by definition of S, M.\ C S. 

Then, one can connect x (resp. x') to a point in Bi V\M.\ (resp. Bj n A4i). 
This shows that one can suppose without loss of generality that x £ JBjD Ai\ 
and x' G nMi. 

Let 7 : [0, 1] — > C be a semi-algebraic path that connects x to x', and let 
G = 7 ^(CnM) and H = [0, 1] \ G. 

Since M\ is finite, we can assume without loss of generality that G is a 
finite set of points, and H is a union of a finite number of open intervals. 

Since 7(G) C M.\ C 5, it suffices to prove that if i and t' are the end 
points of an interval in H, then j(t) and 7(£ / ) are connected by a semi- 
algebraic path inside CDS. 

Notice that j((t,t')) nM\ = 0, so that j(t) and j(t') belong to the same 
Bg by Lemma 13.71 2 b). Recall now that 7(f) and 7(t') both lie in 1?^ n 5 
and that -B^ n S is semi-algebraically connected by (a). Consequently, j(t) 
and "f(t') can be connected by a semi-algebraic path inf^ncScCritS. □ 

We are going to need the following corollary. 

Corollary 3.9. Let 

(Vjp, W, {Mi)x<i<i, (2?i)i<i< 2 ) 

satisfy Property \3. SI 

M = THj,](^i),M = ^l,p](-^2), 

and C R p a finite set containing A/i U A/2. For every semi-algebraically 
connected component C of V, 

c n (w u w) 

is non-empty and semi-algebraically connected. 

Proof. Follows immediately from Proposition 13.31 and Property 13.21 b). □ 

4. Block representations and curve segments 

We denote by D an ordered domain contained in a real closed field R and 
by C the algebraically closed field R[i] . All the polynomials in the input and 
output of our algorithms have coefficients in D and the complexity of our 
algorithms is measured by the number of arithmetic operations (addition, 
multiplication, sign determination) in D. 

In this section, we first define certain representations of points, as well as 
semi-algebraic curves, that are going to be used in the inputs and outputs of 
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our algorithms. Several of these representations share the common property 
that a certain initial number of coordinates are fixed by a triangular system 
of equations, along with certain Thorn encodings and the remaining coordi- 
nates are defined by rational functions to be evaluated at a fixed real root 
of another polynomial (see Definitions 14.11 and 14.81 below). The structure of 
these representations reflect the recursive structure of our main algorithms 
described in Section [71 

After defining these representations we recall the input, output and com- 
plexity of a key algorithm, Algorithm[TJCurve Segments), which is described 
in full detail in [3]. Algorithm Q] accepts as input a polynomial defining a 
bounded real algebraic variety (with some coordinates fixed by a triangular 
system as mentioned above), and outputs a semi-algebraic partition of the 
first (non-fixed) coordinate, as well as descriptions of semi-algebraic curve 
segments (as well as points) parametrized by this coordinate satisfying cer- 
tain properties - which are key to the construction of the main roadmap 
algorithm. Indeed, the curve segments appearing in the main roadmap al- 
gorithm would be limits of the curve segments output by the various calls 
to Algorithm [TJ 

We begin with a few definitions. 

Definition 4.1. A Thorn encoding f,a representing an element a G R 
consists of 

(1) a polynomial / £ D[T] such that a is a root of / in R, 

(2) a sign condition a on the set Der(/) of derivatives of /, such that a 
is the sign condition satisfied by Der(/) at a. 

Distinct roots of / in R correspond to distinct Thorn encodings [21 Propo- 
sition 2.28]. 

A real univariate representation g,r,G representing x G R fc consists of 

(1) a Thorn encoding g, r representing /3 G R, 

(2) G = {go,gi, ■ ■ ■ ,gk) G D[T] fe+1 where g and go are co-prime and such 



4.1. Block representations. In our algorithms, we make recursive calls, 
where we fix blocks of several coordinates. This makes necessary the follow- 
ing rather technical definitions. 

Definition 4.2. A triangular Thorn encoding T = (/m, • • • , f\m])i a repre- 
senting t = (ti, . . . , t m ) in R m consists of 

(1) a triangular system T = (/pj, . . . , / [m] ), i.e. fa G D[Ti, . . . ,Tj] for 
i = 1, . . . , m, such that the zero set of J- in C r is finite; 

(2) a list, a = (ai, . . . , o~ m ), where for i = l,...,m, o~i is the Thorn 
encoding of the root U of /m(ti, . . . , ij-i, Ti). 

A triangular Thorn encoding is quasi-monic if the leading coefficient of /m G 
D[T"i, . . . , Tj] with respect to Tj is a strictly positive element in D. 



that 
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Let P, Q be two polynomials in D[T] with D a domain, the pseudo- 
remainder of P by Q with respect to T is defined as 

PsRem T (P,Q) = Rem T {bP~ q ~ 1 P, Q) 

where q = deg T (Q) and the leading coefficient of Q with respect to T is b q . 
Note that, with 

P = CQ + PsRem T (P,Q) 

both C and PsRemy(P, Q) have coefficients in D. 

The pseudo-reduction modulo F associated to / G D[Ti, . . . , T m ] is defined 

as 

(4.1) PsRed(/,.F) = PsRem Tl (. . . (PsRem Tm (/, 9[m] ), . . .), 9[1] ). 

The pseudo-reduction involves only coefficients in the domain D since the 
triangular Thorn encoding is quasi-monic and pseudo-division is used. Note 
that at the zeros of F the signs of / and PsRed(/, F) coincide. 

Remark 4.3. If / € Dpi, . . . ,T m ] is a polynomial of degree D, and d is a 
bound on the degree of the fa with respect to each T, the complexity of 
computing PsRed(/, F) is (Dd)°^ (see Section [8J Proposition [53])). 

Definition 4.4. A real block representation F, a, L, F representing y € 
consists of 

(1) a triangular Thorn encoding F = (/m, • • • , f\ m ]), & representing a 
root t = (ii, . . . , t m ) of F in R m ; 

(2) a list of natural numbers L = ,£ m ) such that 

£ = h + ---+£ m ; 

(3) a list of polynomials F = (-Fji], ■ ■ ■ , -F] m ]), where 

F[i\ = (/[i]o> • • • » /[i]4)> G D I T i> • • • » ^i], < j < £i, 

with f[i](h, ti- 1} Ti), f[i]o(t 1} Ti) coprime (as polynomi- 
als in Ti), such that 

tt 



y = (y[i],- ■ ■ ,y[ m ]) e R > 



with 



f f{i\i(h, ■ ■ ■ ,U) f[i\£i(ti,---,k)\ 
V[{\ = 7 — t, rr> ■ • • > -7^77 7T > 1 < * < ™- 

V/[i]0(*l, ■ ■ ■ ,U) J[i]o[tl, ■ ■ ■ ,U) J 

In case l\ = • • ■ = £ m = p, then we will write 
(4.2) L=[p m }. 

Notation 4.5 (Substituting a real block representation in a polynomial). 
Let F , a, L, F be a real block representation representing y £ R^, and t £ R m 
represented by J 7 , a. 
Let 

■• • , T i) f[i]£.(Ti, ...,Ti) 



f li] (T 1 ,...,T i )=( 



f\i]o(Ti, ■ ■ ■ ,Ti) ' ' /[j] (Ti, . . . , Ti) 
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Given Q £ D[Xl, . . . ,X k ] with I < k, we define T = (T\, . . . ,T m ), Qp G 
B[T,X £+1 ,...,X k ] by 

(4.3) Qf ■= fo{T)Q (jf[i](?i), . . . , f[ m ](Ti, . . . ,T m ),X i+ i, . . . ,X k ) , 
where 

m 

/o(T) = n/[i]o(ri > ...,T i ) e S 

i=l 

and e, is the smallest even number > deg^ (Q), where X^ is the block of 
variables Xi 1+ ... + £._ 1+ i, Xi 1+ ... + £ i . 
Note that 

Q F (t,Xe+i, . . . ,X k ) = f (t)Q(y,Xi + i, . . . ,X k ), 
with f {t) > 0. 

Notation 4.6 (Substituting a real block representation in a parametrized 
univariate representation). Let J-,o,L,F be as above and let 

g,G = (g ,ge + i, . . .,9k), 

be a parametrized univariate representation with g,gi G D[Xi, . . . ,Xg, U], 
where X±, . . . , Xi are the parameters. 

We denote by G F the tuple (go,F, ■ ■ ■ , 9k,F), where each g^p S D[T, U] 
and is defined by 

(4.4) g it p := f {T)gi (/[i](?i), . . .,/[ m ](Ti, . . . ,T m ),X i+1 , ...,X k ) , 
where 

/o(r) = n/ W o(Ti,...,T i )% 

i=l 

and ej is the smallest number > max.,- deg X[i] (gj) ■ 

Definition 4.7. Let i G R m be represented by a triangular Thorn encoding 
J 7 , a. 

A Thorn encoding g, r representing (3 over t consists of (using the same 
notation as above) 

(1) a polynomial g G D[Ti, . . . , T m , T] such that g(t, /3) = 0, 

(2) a sign condition r on Deixig) such that r is the sign condition sat- 
isfied by the set Der<r(<?(t, T)) at /3. 

A reaZ univariate representation representing x £ K k over t, consists of 

(1) a Thorn encoding g,r representing /3 over t, 

(2) G = (jg ,g 1 ,...,g k ) G D[T 1; . . . , r m , C/] fc+1 such that g(t,U), g (t,U) 
are coprime, and such that 

= ( gi{t,p) g k (t,P) \ Rfc 
U>M)'''''a>M)/ 
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A real univariate representation over t is quasi-monic if the leading monomial 
of g with respect to U is in D. 

A triangular Thorn encoding representing z = (zi, . . . , z r ) over t with 
variables Y\ , . . . , Y r consists of 

(1) a triangular system % = (hi, . . . , h r ), with 

hi G Dpi, . . . , T m , Yi,...,Yi] 

for i = 1, . . . ,r, such that the zero set of H(y, Z\, . . . , Z r ) in C r is 
finite; 

(2) a list, p = (pi, ... , p r ), where for i = 1, . . . , r, pi is the Thorn encod- 
ing of the root z\ of hi(t, zi, ... , Zi-i,Yj). 

4.2. Curve segments. 

Definition 4.8. Let t G R m be represented by a triangular Thorn encoding 
J-, a. A curve segment with parameter Xj overt on (a.1,01.2) in R fc , 

/i,cn> f 2 ,o-2,g,T,G 

is given by 

(1) cti, ct2 € R represented by Thorn encodings fi, a\ and f 2 , 02 over t; 

(2) a parametrized univariate representation with parameter Xj, i.e. 

9,G = (go,gi, ■ ■ ■ ,9k), 

with gj = Xjg and g,g , ...,g k 'm D[T U . . .,T m ,Xj,U]; 

(3) a sign condition r on Der{/(g) such that for every Xj G (ai, 02) there 
exists a real root u(xj) of g(t,xj,U) with Thorn encoding r, and 
g (t,Xj,u(xj)) + 0. 

The curve represented by f±, o~i, f 2 , o~2, g, r, G is the image of the smooth 
injective semi-algebraic function 7 which maps a point x,- of (01,02) to the 
point of R fc defined by 

f x A = ( 9i(t,Xj,u(xj)) g k (t,Xj,u(xj)) \ 
3 \go(t,Xj,u(xj))' '"' " '' g Q (t,Xj,u(xj)) J ' 

Let Q G R[Xi, ...,X k ). For < I < k and y G R f , we denote 

(4.5) Q(y, -) = Q(yi, . . .,y e ,X e+1 , . . .,X k ). 

Remark 4.9. Abusing notation slightly, we will occasionally identify 

ZeT(Q(y,-),R k - £ )cR k - £ 

with 

{y} x Zer(Q(y, -), R fc ^) = Zer(Q, R k ) y C R fc , 

and more generally, for a semi-algebraic set A C R k , A y C R k with {x G 
R k - e \(y,x)GA y }. 

We now recall the input, output and complexity of O Algorithm 15.2 
(Curve segments)]. 
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Algorithm 1. [Curve Segments] 

Input. (1) a point t G R m represented by a triangular Thorn encoding J 7 , a; 

(2) a polynomial P G DfTi . . . , T m , Xi, . . . ,X k ] for which Zer(P(t, -), R 
is bounded; 

(3) a finite set of points contained in Zer(P(i, — ), R fc ) represented 
by real univariate representations U over t. 

Moreover, all the polynomials describing the input are with coef- 
ficients in D. 

Output. (1) An ordered list of points c\ <...< cjv of R with Ci, i = 1, ... ,N 
represented by a Thorn encoding g%,Ti over t. The Cj's are called 
distinguished values. 

(2) For every i = 1, . . . , N, a finite set of real univariate representa- 
tions T>i with parameter X\ over t representing a finite number 
of points, called distinguished points. 

(3) For every i = 1, . . . , N — 1 a finite set of curve segments Ci 
defined on (cj, Cj+i) with parameter Xi, over t. The represented 
curves are called distinguished curves. 

(4) For every i = 1, . . . , N — 1 a list of pairs of elements of Ci 
and T>i (resp. 2?j+i) describing the adjacency relations between 
distinguished curves and distinguished points. 

The distinguished curves and points are contained in Zer(P(i, — ), R fc 
The sets of distinguished values, distinguished curves, and dis- 
tinguished points satisfy the following properties. 
CSi. If c G R is a distinguished value, the set of distinguished 

points output intersect every semi-algebraically connected 

component of Zer(P(t, c, — ), R fc_1 ). 

If c G R is not distinguished, the set of distinguished 
curves output intersect every semi-algebraically connected 
component of Zer(P(t, c, — ),R fc ). 
CS2- For each distinguished curve output over an interval with 
endpoint a given distinguished value, there exists a distin- 
guished point over this distinguished value which belongs 
to the closure of the curve. 
Complexity. If d = deg x (P) > 2, deg T (P) = D, and the degree of the polyno- 
mials in J- and the number of elements of IA are bounded by D, the 
number of arithmetic operations in D is bounded by £) ( m )(l ( mk ) . 
Moreover, the degree in Tj of the polynomials appearing in the out- 
put is bounded by Dd°( k \ 



5. LOW DIMENSIONAL ROADMAP IN A SPECIAL CASE 

In this section we describe an algorithm for computing the roadmap of 
a variety described by equations having a special structure. Although, 
this algorithm is very similar to [SJ Algorithm 15.3 (Bounded Algebraic 
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Roadmap)], the complexity analysis differs because of the special structure 
assumed for the input. 

Let Q G R[Xi, . . . , Xjf] and suppose that V = Zer(Q, R k ) is bounded. 

For < £ < k, < p < k — £, and y G R e , we denote 

dQ(y,-) dQ(y,- 



(5.1) Cr p (Q(y,-)) d = f (q{ V ,-), 



9Xg +p+ 2 dXk 
We assume that Q satisfies the following property. 

Property 5.1. For every £, < £ < k, < p < k — £, and y G K e , the 
algebraic set 

WV = Zer(Cr p (Q(y,-)),R k - e ) 

is of dimension p. 

Remark 5.2. Note that for every y G R^, z G R r , (Wy) z = W9 y ^ has a 
finite number of points and intersects every semi-algebraically connected 
component of V( ViZ \ by [H Proposition 7.4] . 

Now suppose that Q satisfies Property 15.11 For every p, 1 < p < k, and 
y G R^, with £ + p < k we are going to define (My,i)i<i<2i (Py,i)i<i<2 so 
that the tuple 

(V y ,p, W$, {My,i)l<i< 2 , (p y ,i)i<i<2) 
satisfies Property 13.21 using a slight abuse of notation (cf. Remark I4.9p . 

Definition 5.3. Let 

(1) M y ,i = W° C Vy, and V V)1 = vr m (VF°); 

(2) Vy t 2 C R the set pseudo-critical values (see [21 Definition 12.41]) of 
-Kt+i on 7T[£ +1) fc](Wy ) and M. Vl 2 a set of points such that for every 
c G T>2, M2 intersects every semi-algebraically connected component 
D of (W^) c . 

Note that by Property 15.11 Wy is of dimension p (in particular, W® is fi- 
nite), and satisfies Propertv l3.2l 2): for each z G R p , (Wy) z = W? z ^ is a finite 
set of points having non-empty intersection with every semi-algebraically 
connected component of V( VjZ \ by Remark [5.2i Moreover, Zer(Q(y, — ),~R k ~ i ) 
is clearly bounded (since Zer(Q, R fc ) is bounded), and the finite set M y> i = 
W° is the union of the Q(y, — )-singular points of Zev(Q(y, — ), R fc and the 
Q(y, — )-critical points of the map tt£ + i on V y = Zer(Q(y,-),R fc ^). Thus, 
M y ,i and V Vt i satisfies Property 13.21 1). 

Note also that M y ,2 satisfies Property 13.21 3). Indeed, the intersection of 
M.y,2 with every semi-algebraically connected component of (Wy) c is non- 
empty, by [21 Proposition 12.42]. Moreover for every interval [a, b] and c G 
[a, b] such that [a, b] contains no point of T> y> 2, except maybe c, and for 
every semi-algebraically connected component D of Wi y y x \ a u, D( y ^ c \ is a 
semi-algebraically connected component of (W y ) c , by [21 Proposition 15.4]. 

So we have proved 
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Proposition 5.4. If Q satisfies Property 15, 11 using the notation above, 

(Vy,P, Wjj, (My,i)l<i< 2 , (T>y,i)l<i< 2 ) 

satisfies Property \3. 21 

We are going to describe below, in the special case where Q satisfies 
Property 15.11 an algorithm directly adapted from [2J Algorithm 15.3] for 
computing a roadmap of certain subvarieties of Zer(Q,R fc ) of dimension 
at most p : this is Algorithm [2] (Roadmap for Lower Dimensional Special 
Algebraic Sets) 

Remark 5.5. In all our algorithms, the roadmaps output are represented by 
a finite number of real univariate representations and curve segments over a 
point defined by a triangular Thorn encoding (see Definitions 14.41 14.71 and 
~8\ above). 



Algorithm 2. [Roadmap for Lower Dimensional Special Algebraic 
Sets] 

Input. (1) a polynomial Q £ D[Xl, . . . , X^] satisfying Property 15. 1\ and 
for which V = Zer(Q,R fc ) C B fc (0, 1/c) (where c G R); 

(2) numbers p, m, r > satisfying mp < k, < r < p; 

(3) y E R mp represented by a real block representation J 7 , a, \p m ],F 
(see (|4.2p ) with t G R m represented by a quasi-monic triangular 
system J 7 , a; 

(4) z£R r represented by a triangular Thorn encoding H,p over t, 
with variables X mp+ i, . . . , Ar m p +r .; 

(5) a finite set of points Mo contained in W?~£\ represented by real 
univariate representations Uo, over (t,z) (using the notation of 
Property I5.ip . 

Output, a roadmap BM(Wf~^, Mo) for (W?~^,Mo) represented as a union 
of of curve segment representations and real univariate representa- 
tions over points defined by triangular Thorn encodings. The adja- 
cencies between the images of the associated curves and points are 
also output. 

Complexity. flOKpJ^tKp)* w here d = deg(Q) > 2, and D is a bound on the 
degree of H, J 7 , F and the number and degrees of the elements in Uo- 

Procedure. 

Step 0. Define 

P := ^2 ^ e D[Ti , . . . , T m , X mp+ \ , . . . , Xk] 

A£Cr p (Q F ) 

using Notation 14.51 and (15. 1|) . and initialize i := 1. 
Step 1. If i — 1 + i < p, call Algorithm [T] (Curve Segments) with input 

(■F,n),(a,p),P,Uo- 
Pseudo-reduce modulo T using (|4.1h and place the output in the 
description of RM(W?~J,.Mo). 
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Step 2. Set i := i + 1. Using the notation in the output of Algorithm [H for 
every j = 1, . . . , N , define 



z : = 


(Z,Cj), 




n:={n, gj {T x ,... 






p ■= 






U : 


= v 3 , 





and call Step 1 of Algorithm [2] (Roadmap for Lower Dimensional 
Special Algebraic Sets) recursively, with input 



Proof of correctness. Notice that 

Wf-Jj = Zer(P(i,s,-)),RM"*4-r)). 

The correctness of the algorithm then follows from the correctness of Al- 
gorithm Q] (Curve Segments) and Proposition 12.21 The only additional fact 
that needs to be checked is that when the recursion ends with r = p, the 
algebraic variety Zer(P((i, z, z'), — ), R fc_ p( m + 1 )) is zero-dimensional, where 
z' = (cj 1 , . . . ,Cj p _ r ) S R p_r and the various Cj i G R are associated to the 
Thorn encodings computed in Step 1 of the algorithm. This is the case 
because Zer(P((t, z, z'), -), R k ~p( m + 1 )) = W? z ,y and W? zl) is zero- 
dimensional by Property 15.11 □ 

Complexity analysis. The depth i of the recursion is bounded by p — r, 
and the total number of recursive calls at depth i is bounded by d ^ lk \ 
Thus, there are at most d 0<y<yP ~ T ^ calls to Algorithm [1] (Curve Segments). 

In each of the calls to Algorithm [1] (Curve Segments), the number of 
arithmetic operations in D is bounded by L)°( m +P)(i c, (( m +P) fc ) using the com- 
plexity analysis of Algorithm [1] ( Curve Segments). Moreover the number of 
arithmetic operations needed for each pseudo-reduction is (Dd k )°^ m ^ since 
the degree in Tj of the output of Algorithm [T] (Curve Segments) is DdP^ 
using Remark 14.31 

Thus, the total number of arithmetic operations in D for Algorithm [2] is 
bounded by £)0(m+p)^0((m+ P )fe)_ n 



6. Low DIMENSIONAL roadmap in general 

In this section, we first explain how to perform an infinitesimal defor- 
mation of any given polynomial Q 6 RLXi, . . . , X^] such that the deformed 
polynomial satisfies Property 15.11 

We then sketch how to compute the limit of a curve, and finally how to 
compute the limits of roadmaps of certain algebraic sets which are the critical 
locus of dimension p of certain projection maps restricted to the algebraic 
hypersurfaces obtained after performing an infinitesimal deformation. 
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6.1. Deformation. We consider a bounded algebraic set denned by a non- 
negative polynomial Q. Our aim is to define a deformation of Q defining a 
polynomial satisfying Property 15.11 

Suppose that the polynomial Q G R[Xl, . . . , Xk], and the tuple {d\, . . . , dk) 
satisfy the following conditions: 

(1) Q(x) > for every x G R k , 

(2) Zer(Q, R k ) is bounded, 

(3) di > d 2 ■ ■ ■ > d k , 

(4) deg(Q) < d x , tDeg x .(Q) < d u for i = 2, . . . , k. 

Let di be an even number > di, i = 1, . . . , k, and d = [d\, . . . , dk)- Let 

G k (J) = 1 + • • • + + X 2 2 + • • • + X 2 k + + 2k, 

and note that V x Gk{d){x) > 0. 
We denote 

Notation 6.1. 

(6.1) Def(Q,e) = -eG k {d) + Q, 

(6.2) V y>£ = Zer(Def(Q,e)(y,-),R(e) k+1 ) 

/f% q\ r (n v ^, rn v aDef(Q,e) dDef(Q,e) 

(6.3) Cr^(Q,e) = Def(Q,e), 



(6.4) M", = Zer(Cr p (Q,e)(y,-),R<e 



A;+l^ 



Proposition 6.2. For every t, < i < k, and every y £ R , 

a) Def (Q, e)(y, — ) satisfies Property \5. 11 

b) lim e induces a 1-1 correspondence between the bounded semi- algebraically 
connected components of 

V y , e = Zer(Det(Q,e)(y,-),R(e) k+1 ) 

and the semi- algebraically connected components of 

Z y = Zer(Q(y,-),R k ) 

Proof, a) follows from [21 Proposition 12.44] and b) from [2l Lemma 15.6]. 

□ 

We are going to describe in Section 16.31 an algorithm for computing the 
limit, under the lim e map, of a roadmap of the critical locus of dimension p, 
Wy j£ , of Vy t£ = Zer(Def(Q, s){y, —),R(e) k+1 ~ £ )- In order to achieve this we 
first need to compute limits of curves, which is the purpose of Section | 



6.2. Limits of points and curve segments. The general problem of com- 
puting the image of a semi-algebraic set 5 C R(e) k which is bounded over 
R under the lim e map reduces to the problem of computing the closure of a 
one-parameter family of semi-algebraic sets, which can be done using quan- 
tifier elimination algorithms (see, for example, [21 pg. 556]). However, the 
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complexity of this general algorithm, d k , is not good enough for our pur- 
poses in this paper. Fortunately, we need efficient algorithms for computing 
limits only in two very special situations, where we can do better than in 
the general case. 

These two special cases are the following : 

(1) when the set is a point represented by a real univariate representa- 
tion, 

(2) when the set is a curve represented by curve segments. 

We give now the input output and complexity of Algorithm [3] (Limit of a 
Bounded Point) and Algorithm H] (Limit of a Curve). A full description of 
these algorithms, their correctness and complexity analysis appear in Section 

El 

Algorithm 3. [Limit of a Bounded Point] 

Input. (1) a quasi- monic triangular Thorn encoding J 7 , a, with coefficients 
in D, representing a point t G R m ; 
(2) a real univariate representation g e ,T e , G £ over t with coefficients 
in D[e], representing a point z £ G R(e) p bounded over R. 
Output. (1) a quasi-monic triangular Thorn encoding J-', a', representing 
the point t G R m ; 
(2) a quasi-monic real univariate representation (h, H) representing 

z = lim z e € R p . 

e 

Complexity. If D\ (resp. D2) is a bound on the degrees of the polynomials in 

J 7 , g £ and G £ with respect to T±, . . . , T m (resp. e, U), then D\ (resp. 

D2) is a bound on the degrees of the polynomials appearing in the 

output, and the number of arithmetic operations in D is bounded by 
D 0(m) D 0(l) 

Remark 6.3. Note that there is a possibly new representation of t in the 
output of Algorithm [3] (Limit of a Bounded Point). The reason for this 
peculiarity is explained in Section [8] (cf Proposition I8.4j) . 

Algorithm 4. [Limit of a Curve] 

Input. (1) a quasi-monic triangular Thorn encoding J 7 , a with coefficients 
in D representing t G R m ; 

(2) a triangular Thorn encoding H £ ,p £ over t with coefficients in 
D[e] representing z £ G R(e) r over t; 

(3) a curve segment with parameter X r+ \ and coefficients in D[e] 
over (t,z £ ), representing a curve bounded over R. 

Output. (1) a real univariate representation p z ,p z ,Pz of z = lim £ (z £ ), with 
u the root of p z with Thorn encoding p z ; 
(2) a finite set {di, . . . , djv-i} where each d{ is a real univariate 
representation over (t, u, q), and q is given by a Thorn encoding 
over t fixing A m(i) ; 
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(3) a finite set {w±, . . . , wjy}, of curve segments over (t, u) with Wi 
parametrized by X^y 

Moreover, the union of the curves represented by W, and the 
points represented by V define a partition of 5 = lim e (S e ). All 
the coefficients of the polynomials in the output belong to D. 
Complexity If the polynomials occurring in the input have degrees bounded by D, 
then the complexity of the algorithm is bounded by D°( m + r ) , 

6.3. Low dimensional roadmap algorithm. We are going to describe 
an algorithm computing the limit of a roadmap of the critical locus of di- 
mension p, Wy j£ , of the deformation V V)£ = Zer(Def (Q, e)(y, — ), R(e) k+1 ~ e ) 
of Z y = Zer(Q(y, — ), K k ~ e ). The algorithm proceeds by first calling Algo- 
rithm [2] (Roadmap for Lower Dimensional Special Algebraic Sets) in order 
to compute a roadmap for Wy jE , and then computes the image of the result- 
ing roadmap under the lim e map. Note that this limit is not necessarily a 
roadmap of V y , since a semi-algebraically connected component of V y might 
contain the image under lim e of more than one semi-algebraically connected 
components of Wy >£ . 

Algorithm 5. [Limit of Roadmaps of Special Low Dimensional Va- 
rieties] 

Input. (1) a natural number p < k; 

(2) a polynomial Q G D[Xi, . . . ,X k ] for which Z = Zer(Q,R fc ) C 
B fc (0,l/c) (with c G R); 

(3) y G R mp represented by the real block representation 

(see (|4.2p ) with coefficients in D, such that t G R m is represented 
by a quasi-monic triangular Thorn encoding J-, a; 

(4) a finite set of points Afy t£ G {y} x R(e) p represented by quasi- 
monic real univariate representations V e , over t. 

Output. Real univariate representations and curve segment representing the 
set of points 

K = (7r [1)fe] o hm)(RM(^ e , W My J) 
where Wy >e is the zero set of d p {Q{y, — ),e), 



(6-5) W My , e = {Wl e ) Mv , e = |J W, 



o 



and RM« £ ,% W ) is a roadmap for (W^ £ ,Wm v , e ) 
Complexity. j DO(m+ P ) d o((m+ P )fc where d = d eg (Q) > 2 and D is a bound on the 
degree of J 7 , F and the number and degrees (including that in e) of 



Procedure. 



the elements in M V)£ . 
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Step 1. Let T = (T\, . . . ,T m ), and, using (14. 5p and (I6.3D . 

P= J2 A 2 eD[e,T,X mp+1 ,...,X k ]. 

A£Ct p (Q f ,s) 

Call [21 Algorithm 12.18 (Parametrized Bounded Algebraic Sam- 
pling)] with input P and parameters e,T, X mp+ i, . . . ,X( m+1 \ p and 
output a set of parametrized univariate representations with variable 
U. 

Pseudo-reduce them modulo T using (|4.ip and place the result in 
U>. 

For every (h e ,H e ) S U e , and every z £ £ Ny,s represented by a 
real univariate representation (g e ,r,G E ) £ V e , use [21 Algorithm 
12.20 (Triangular Thorn Encoding)] with input the triangular sys- 
tem {F,g e ,h e ) to compute the Thorn encodings of the real roots of 
h e (y,z £ ,U). Let W y Ze be the set of real univariate representations 
over y, z e so obtained. Define 

U 'y,e = (J U 'v,ze- 

The set of points represented by U' y e is Wj\f y e (see (16. 5p ). 

Step 2. Call Algorithm [2] (Roadmap for Lower Dimensional Special Alge- 
braic Sets) with input Def(Q,e) (see (|6.1j) ) and p, the real block 
representation J 7 , a, [p m ],F, r := and W y E . The output of Algo- 
rithm [2] (Roadmap for Lower Dimensional Special Algebraic Sets) 
consists of a set of real univariate representations and curve seg- 
ments over triangular Thorn encodings. Each such curve segment, 
7 = (/i, C7i, /2, 0"2j 9i t, G), is defined over some (t,z-y) with r 7 < p 
and z 7 G R(e) r7 , represented by a triangular system F^Tl^. 

Step 3. For each such curve segment 7 over (t, Zj), output in the previous 
step over call Algorithm [5] (Limit of a Curve) with input the trian- 
gular system J 7 , % 7 and 7. Finally, project on R fc by forgetting the 
last coordinate. 

Remark 6.4. The role played by the set of points Wj\f E which are included 
in the roadmap of Wy }£ , whose limit is computed by Algorithm [5] (Limit 
of Roadmaps of Special Low Dimensional Varieties), will become clear in 
the proof of correctness of Algorithm [6] (Baby-giant Roadmap for Bounded 
Algebraic Sets) (see £L2}). 

Proof of correctness. First note that it follows from Proposition 16.21 
that Def(Qi?,e) satisfies Property 15.11 and hence (Wy >e )sf e is a finite set of 
points. The correctness of the algorithm now follows from the correctness of 
Algorithm [2] (Roadmap for Lower Dimensional Special Algebraic Sets) and 
Algorithm [4] (Limit of Curve) . □ 
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Complexity analysis. The number of arithmetic operations performed 
in D[e] in Step 1 is bounded by D°( m +P) d°(( m +p) k ) arithmetic operations in 
D[e] according to the complexity analysis of [21 Algorithm 12.18 (Parametrized 
Bounded Algebraic Sampling)] and [21 Algorithm 12.20 (Triangular Thorn 
Encoding)]. Since the degree in e in the output of [21 Algorithm 12.18 
(Parametrized Bounded Algebraic Sampling)] is and does not change 

during the pseudo-reduction, the number of arithemtic operations in D in 
Step 1 (and hence the complexity) is bounded by D°( m +p)d°(( m +p) fc ). 

The number of arithmetic operations performed in D[e] in Step 2 is 
bounded by D°( m +!')(i C( m +P) fc ) according to the complexity analysis of Al- 
gorithm [2] (Roadmap for Lower Dimensional Special Algebraic Sets). More- 
over the degree in e is bounded by 0(d) k by [21 Algorithm 15.10 (Parametrized 
Curve Segments)], since the computations of Algorithm (Curve segments) 
with coefficients in D[e] is contained in the computations of [21 Algorithm 
15.10 (Parametrized Curve Segments)] in D with parameter e. So the the 
number of arithmetic operations in D in Step 2 (and hence the complexity) 
is bounded by L>0(m+ P ) d o((m+p)fc)_ 

The complexity of Step 3 is also bounded by D°( m +p) d°(( m +p) k according 
to the complexity analysis of Algorithm U] (Limit of Curve). 

Thus the total complexity of the algorithm is D°( m+ ^ d oi -( m+p ^ . n 



7. Main result 

We now describe our main result Algorithm [7J (Baby-giant Roadmap for 
General Algebraic Sets). It is based on Algorithm [6] (Baby-giant Roadmap 
for Bounded Algebraic Sets), computing a baby step - giant step roadmap al- 
gorithm for a bounded algebraic set. The algorithm for computing roadmaps 
of general (i.e. not necessarily bounded) algebraic sets, Algorithm [7] (Baby- 
giant Roadmap for General Algebraic Sets) is then obtained from Algorithm 
[6] (Baby-giant Roadmap for Bounded Algebraic Sets) following a method 
similar to the one in [2] to go from the bounded case to the general case. 

Algorithm [6] (Baby-giant Roadmap for Bounded Algebraic Sets) proceeds 
roughly as follows. We denote by y the fixed coordinates. If the number 
of non-fixed coordinates is too small (i.e. less than the number p which is 
prescribed in the input), then we compute the roadmap using [21 Algorithm 
15.3 (Bounded Algebraic Roadmap)]. Otherwise, we compute representa- 
tions of points in J\f y C W defining the fibers at which we make recursive 
calls to the same algorithm; these are the giant steps. 

For the baby steps, the algorithm uses Algorithm [5] (Limit of Roadmaps 
of Special Low Dimensional Varieties) to compute the limit (under the lim e 
map) of the roadmap of the critical set Wy )£ going through a well chosen 
finite set of points. 

We are now ready to proceed to the description of Algorithm [6] (Baby- 
giant Roadmap for Bounded Algebraic Sets) 
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Let as in Notation 16,11 

V y , £ = Zer(Def(Q,e)(y,-),R(e) fc+1 ) 
W£ £ = Zev(Cv p (Q,e)(y,-),R(e) k+1 ) 

and define 

{■M-y l£l i)l<i<2, \Dy,e,i)\<i<2) 

from Vy t£ ,Wy }£ as in Definition 15.31 

It follows from Proposition 16.21 and Proposition 15.41 that 

{Vy, S ,P,Wy >e , (My, £ ,i)l<i<2, (£W)i<i<2) 
satisfies Property 13. 21 

Algorithm 6. [Baby-giant Roadmap for Bounded Algebraic Sets] 

Input. (1) a polynomial Q G D[X l7 . . . ,X k ] which Z = Zer(Q,R fc ) C 
B fe (0,l/c) (where c G R); 

(2) y G R mp represented by a real block representation 

(see (|4.2p ) such that t G R m is represented by a quasi-monic 
triangular Thorn encoding F, a; 

(3) a finite set of points M y fl in Z y = Zei(Q(y, —),R k ~ mp ) repre- 
sented by quasi-monic real univariate representations Uq over t. 
All the coefficients of the input polynomials are in D. 

Output, a roadmap, BGRM(Zy, M y fi), for (Z y ,M y fl)- 
Complexity. d°( k2 /p+P k ) operations in D where d = deg(Q) > 2 and the degrees 
of the polynomials in J 7 , F, as well as the degrees of the polynomials 
and the number of elements in IAq are all bounded by d°^ k \ 

Procedure. 

Step 1. If (m+l)p > k call [2J Algorithm 15.3 (Bounded Algebraic Roadmap)] 
with input 

(1) the quasi-monic triangular Thorn encoding J 7 , a representing 
t G R m , 

(2) the polynomial Qf, using Notation 14.51 

(3) the finite set of points M y ,o in Z y = Zev{Q F {t, -),R k ~ mp ) rep- 
resented by real univariate representations Uq over t. 

Otherwise set i := 1 and do the following. 
Step 2. Determination of the finite set of points M y used in the recursive 
call. Let T = (Ti, . . . , T m ), and, using (j4.5j) and (|6.3p . 

P= A 2 eB[e,T,X mp+1 ,...,X k }. 

AeCr p (Q F ,e) 

Step 2 a). Call [21 Algorithm 12.18 (Parametrized Bounded Algebraic Sam- 
pling)] with input P and parameters e, T, and output a set of parametrized 
univariate representations with variable U. Pseudo-reduce them 
modulo F using (14. lh and place the result in U e ^\. 
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For every [h £ ,H £ ) G U £ ,i-, use [21 Algorithm 12.20 (Triangular 
Thorn Encoding)] with input the triangular system (J 7 , h £ ) to com- 
pute the Thorn encodings of the real roots of h e (y,U). Let Uy )£i i 
be the set of real univariate representations over y so obtained. Let 
M.y j£ ,i CC Vy i£ be the set of points represented by U y , £j i. 

Projecting Uy j£; i, by forgetting its last components, obtain a set 
of quasi-monic real univariate representations V y<e ,i representing 

over t. Then apply Algorithm [3] (Limit of a Bounded Point) with 
Vy je ps,i as input to obtain a set of quasi-monic real univariate rep- 
resentations V Vt i representing 

N y ,i = lim(A^ y , e ,i) 

over t. 

Step 2 b). Perform AlgorithmQ] (Curve Segment) with input P and the triangu- 
lar Thorn encoding J-, a and retain the set of univariate representa- 
tions, Uy^^i representing A4 y>£j 2 C V y>£ , which are the distinguished 
points in the output. 

Projecting Wy, e ,2j by forgetting its last components, obtain a set 
of real univariate representations Vy <e ,2 representing 

•A/y, e ,2 = 7r[mp+l,(m+l)p](^W,e,2)- 

Then apply Algorithm [3] (Limit of a Bounded Point) with Vy j£ 2 as 
input to obtain a set of quasi-monic real univariate representations 
V y ,2 representing 

Af Vt 2 = lim(A/" a)£) 2). 

Step 2 c). Projecting Uq, by forgetting its last components, obtain a set of 
quasi-monic real univariate representations V y $ representing 

■A/j/,0 = 7r [mj)+l,(m+l)p](-Wj,o) 

over t. 
Let 

My = Ny,0l)Ny,\ UA/^2, 
Vy = V,,,0UV,,,lUV,,,2 J 

A/" a , e = J\f y U U A/"y, £ ,2, 

and 

V y , e = V y U V„, e ,i U V„, e , 2 . 
Step 3. Call Algorithm [5] (Limit of Roadmaps of Special Low Dimensional 
Varieties) with input p, Q, the real block representation J 7 , a, \p m ],F, 
and V Vj£ and note that it contains 

WVv,. = = U W Le- 
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Recursive call For every element u = ((J 7 , h), (<t,t), (F, H)) G V y , representing 
(y, z) £ R( m+1 )P ; compute a set, W^z, of quasi-monic univariate rep- 
resentations over ((J 7 , h), (a, r), (F, iJ)), representing 

(7-1) M^z) = (n [1M o\im)W° Zj£ ). 

using Algorithm [3] (Limit of a Bounded Point). 
Call Algorithm [6] (Baby-giant Roadmap for Bounded Algebraic 
Sets) recursively with input 

Q, F := OF, h),a:= (a, r), L := F := {F,H),i :=i + l 

and U^y^)^ '■= M(y,z) U (My,o)zi where (Z^o)z is a set of quasi-monic 
real univariate representations representing (Ai y fl)z- 

Output the set of curve segments computed in the last two steps. 

Remark 7. 1 . Algorithm [6] would have been much simpler if we could make 
recursive calls to Algorithm [6] at the fibers over the points in A/^ j£ , and 
thus obtain a roadmap first of Vy >£ , and finally take the image of the re- 
sulting roadmap under the lim £ map. In this case the proof of correct- 
ness of the algorithm would be an immediate consequence of the main 
connectivity result, Corollary 15.41 and the fact that the image under lim e 
of a bounded, semi-algebraically connected semi-algebraic set is also semi- 
algebraically connected. 

However we are unable to compute limits of semi-algebraic curves given by 
curve segments over a real block representations depending on e with number 
as well sizes of the blocks bounded by 0(y/k) with complexity (j, ( k ^)^ 
because we would obtain a degree m e. 

We overcome this difficulty by making recursive calls to Algorithm [6J not 
at the fibers over the points in J\f y>e , but at the fibers over M y = lim e (jV ei3/ ), 
so that the algebraic sets specified in the input to the various recursive calls 
are then Z^ y z ^ for z € N y . In this approach, the only limits of curve segments 
that are computed are those of the roadmap of W y>e , and we can compute the 
limits of these curve segments without spoiling the complexity, as they are 
not defined over real block representations depending on e. However, since 
the recursive calls are made with fibers of Z y (instead of V y , e ), Corollary 15.41 
is not directly applicable, and we need to be more careful about choosing 
the set of points in the input to the recursive calls. It also makes the proof 
of correctness more complicated. 

Proof of correctness. 
Base case. 

If \(k — mp)/p\ = 1 then the correctness of the algorithm is a consequence 
of the correctness of [21 Algorithm 15.3 (Bounded Algebraic Roadmap)]. 
General case. 

Suppose that \(k — mp)/p] > 1. 
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Denote by BGRM(Z y , A^o) the union of the curve segments output by 
Algorithm [6] (Baby-giant Roadmap for Bounded Algebraic Sets) . 
We have that 

BGRM(Z y ,My fi ) =K y U (J BGRM^^.M^U^o),), 
with 

Ky = (7r [1)fc] o lim)(RM(W* e , W A f V:S )), 
denoting by W^ e the zero set of Cvp(Qp,e) and 

WM v , e = {Wl e )M y , e = |J W^ E , 

(see ([631) and (163]) ). 

Proof of My fi C BGRM(Z J/ , M y>0 ). 

The proof is by induction on \(k — mp) /p\ . 

We suppose by induction hypothesis that for every (y, z) G M y that 

M (y>z) U (M y ,o) z ) C BGBM(Z M ,M iyiZ) U (AV)*)- 
Since, Aiyfl C (Af^cO.z, and by induction hypothesis we have that 

{M yfi ) z ) C BGRM(Z M ,M M U (At^o)^), it is clear that BGRM(Z y ,M yfi ) 
contains -Myfi- 
Proof of RMi. 

The property RMi of BGRM(Z J/ , Aiyfi) is also proved by induction on 
\{k — mp)/p] . 

Let C be a semi- algebraically connected component of Z y , and -D = 
BGRM(Z J/ , -M^) U (Af^o)^) n C. We want to prove that D is semi- 
algebraically connected. 

Suppose that x, x' G D, we are going to prove that there exists a semi- 
algebraic path 7 : [0, 1] — > D with 7(0) = £,7(1) = x' . 
Without loss of generality we can suppose that x G 1Z y . 

Since 

BGRM(Z y ,M (ytZ) U(M y , ) z )=n y U |J BGRM(Z (y>z) ,M { y jZ) U(M y ,o) z ), 
we have that x (resp. x') either belongs to 

Uy 

or to some 

BGBM(Z Ml M M U(A1 V) o),) 

with (y, z) G A/ r 

If x G BGRM^^), M( y ,z) U (-My,o)z) we show that x can be connected 
to a point in A^v) inside 

BGRM(Z M ,^ M U(M, i0 ),) 
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by a semi-algebraic path. It follows from Proposition 15.41 that Def (Q,e)f 
satisfies Property 13.21 (2). and hence we have that W? y z \ £ meets every semi- 
algebraically connected component of Vr VjZ \ e . By [31 Lemma 15.6] each 
semi-algebraically connected component of Z( ViZ \ is image under 7Tm w o lim e 
of a unique semi-algebraically connected component of VL Z ) 6 . It follows 
that each semi-algebraically connected component of Zr ViZ \ meets 

since 

Ka* = TO)* 

Finally, applying the induction hypothesis to BGRM(Z( y 2 ), A / f( J/)2 )U(A / fj /! o)z) 
we have that the intersection of BGRM^^ z ) , M.( y ^ z ) U (My t o) z ) with each 
semi-algebraically connected component of Z( y>z \ is non-empty and semi- 
algebraically connected, and meets M(y tZ y Thus, there exists a semi-algebraic 
path with image in ~BGRM(Z( y z y M-! y , z ) U (-M.y,o)z) joining x to a point in 
M (y,z)- 

Since ■M-( y , z ) C lZ y we can assume that x (and similarly x') is contained 
in TZ y . 

Connectivity when x and x' are contained in 7^. 

Since 

K y = (7r [1)fc] o lim)(RM(W£ £ , WfyJ), 

there exists x £ G RM(W£ e , WjV H , e ) (resp. 4 £ RM(W£ e , WV„ .J) such that 
lim e (x e ) = x (resp. lim e (x' £ ) = x'). 
Let 

5 e = w2; e u(^ je )^, e , 

and C £ the unique semi-algebraically connected component of V VtE such that 
(7T[l,fe] olim e )(C £ ) = C. 

By Corollarv l3.9l since N y ,e,i UA/"y i£i 2 C JV^, e , 5 e nC e is semi-algebraically 
connected. So there exists a semi-algebraic path / y £ : [0, 1] ->5 e nC £ , with 
7e(0) = x e ,7 e (l) = x £ . Moreover, there exists a partition of (0,1) C R(e) 
into a finite number of open intervals and points, such that for every open 
interval I in the partition one of the following holds : 

Case 1: 

c W* te . 

Case 2: there exists z e G Ny,e such that 

7e(J) C V (2/) ^) !£ . 

Since Wy i£ C Vy, e , for each point a G (0, 1) defining the partition 
(7-2) 7e (o) G WV«, 8 C RM(^ )£ ,WV 9 ,J. 

Hence, by definition of M( y>z ) (see (|7.ip ) 
(7.3) 7T[i |fc] o lim)(7 e (a)) G -M^), 
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where lim e (z e ) = z. 

In Case 1, we can replace 7 e (I) by a semi-algebraic path having the same 
endpoints and whose image is contained in 

RM(W£ £ ,WaU 

using RMi for 

RM(W? £ ,W Xy J 

as well as (I7.2p Taking the image under ^[i^] bm e of this new path we 
obtain a semi-algebraic path 

7 : lim(J) — > 7Z y . 

In Case 2, (7r[ l fe ] o lim £ )(7 £ (a)), (7Tp. jfc ] l™e)(7e(&)) both belong to 

BGRM(Z M ,^ M U(M, i0 ) z ) 

using (|7.3p . Using the induction hypothesis for BGRM(Z(j /i2 ), -M^) U 
(My o)z) we have that there exists a semi-algebraic path 

7 : pim(o),lim(6)] -> BGRM(Z (y>2) , M M U (M V)0 )z). 

Finally, we have constructed a semi-algebraic path 7 : [0, 1] — ?■ D with 
7(0) =x, 7 (l) =x'- 

This proves that BGRM^^, -M(y,z) U (-^3/,o)z) H C is non-empty and 
semi-algebraically connected proving RMi. 
Proof of RM 2 . 

Let c G R such that -Z^c) is not empty, and let C be a semi-algebraically 
connected component of Z^ yc y We prove that 

BGRM(Z (v) , U (M y , ) z ) n C 

is not empty. It follows from [3j Lemma 15.6 ] that there exists a semi- 
algebraically connected component, C £ , of VL c ) j£ such that 

C = lirn)(C e ). 

Since C £ is non-empty, let x £ G C £ and let z £ = 7Tr mp+ i ( m+ i) p i(a; e ). It follows 
from Proposition l5.4l that (Wy je ) Ze = W? z \ £ meets every semi-algebraically 
connected component of V( y ^\ >£ . Since (7 £ contains a semi-algebraically 
connected component of V( y>z ^ >£ , we have that 

and thus C £ contains a semi-algebraically connected component of (W y>£ ) c 
(since W y>£ C V y>£ ). Now, since the roadmap 

RM(W£ £ ,WaU 

satisfies RM2, RM(W| > )£ , W/v s e ) has a non-empty intersection with every 
semi-algebraically connected component of (Wy, e ) c , and in particular with 
the one contained in C £ . Taking the image under 7Tm w o lim £ map we get 
that TZ y = (7r[i ! j.]olim £ )(RM(Wy )£ , W^ y e )) has a non-empty intersection with 
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(7r[i )fc ]olim E )((7 E ) = C. Since, BGRM(Z ( j / ^, M(y jZ ) U (M y< o) z ) c contains TZ y , 
this finishes the proof. □ 

Complexity analysis. We first bound the number of arithmetic opera- 
tions in Step 1. Since we assume that the degrees of the polynomials in J 7 , F 
are bounded by d°^ k \ it follows from the complexity analysis of [21 Algo- 
rithm 15.3 (Bounded Algebraic Roadmap)], and O Algorithm 15.2 (Curve 
Segments)], that the number of arithmetic operations in this step is bounded 
by 

d O(km) d O((k-mp) 2 ) _ ^0{km+p 2 ) 

since k — mp < p. 

The number of arithmetic operations in D[e] in in Step 2 is bounded by 

d O(mk) and the de 

gree and number of univariate representations produced is 
bounded by 0(d) k ~ mp . Moreover the degree in e is bounded by 0{d) k by [21 
Algorithm 15.10 (Parametrized Curve Segments)], since the computations 
of Algorithm (Curve segments) with coefficients in D[e] is contained in the 
computations of [21 Algorithm 15.10 (Parametrized Curve Segments)] in D 
with parameter e. So the number of arithmetic operations in D in in Step 2 
is bounded by d°( mfc ). 

The complexity of computing lZ y in Step 3 is bounded by <i (( m +P) fc ) 
given that the number of arithmetic operations of Algorithm [5] (Limit of 
Roadmaps of Special Low Dimensional Varieties) is d°^ m+p ^ k \ 

The total number of recursive calls at depth i is d ^ k% \ and for each 
such call the number of arithmetic operations in D in Steps 1, 2 and 3 is 
bounded by d°(( m + l +p) k +P f where < i < [k/p\ — m. Since the depth 
of the recursion is at most [k/p\ — m, we obtain that the total number of 
arithmetic operations in the domain D is bounded by 

d O(k 2 / P ) d O((k/p+p)k+p 2 ) = d O(k 2 /p+pk)_ 

□ 

We now describe Algorithm!?] (Baby-giant Roadmap for General Algebraic 
Sets) for computing a roadmap of a general (i.e. not necessarily bounded 
algebraic set). This algorithm is essentially the same algorithm as [21 Algo- 
rithm 15.5 (Algebraic Roadmap)], except that we call Algorithm [6] (Baby- 
giant Roadmap for Bounded Algebraic Sets) after reducing to the bounded 
case instead of of [21 Algorithm 15.3 (Bounded Algebraic Roadmap)]. We 
first need a notation. Let P £ Rf-X"] be given by 

P = a p X p H \-a g X q ,p> q,a q a p ^0. 

Notation 7.2. We denote 




c(P) = 

Algorithm 7. [Baby-giant Roadmap for General Algebraic Sets] 
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Input. (1) a polynomial Q G D[Xi, . . . ,X k }\ 

(2) a finite set of points M y fl in Zer(Q,R k ), represented by real 
univariate representations Uq. 
Output, a roadmap, BGRM(Zer(Q, R k ),Mo), for (Zer(Q, R fc ), M))- 
Complexity. ^o( fc2 /p+p fc ) operations in D. 
Procedure. 

Step 1. Introduce new variables X k +\ and e and replace Q by the polynomial 

Q £ = Q 2 + {e 2 {Xl + ... + Xl+ 1 )-l) 2 . 

Replace C R fc by the set of real univariate representations rep- 
resenting the elements of Zer(e 2 (Xf + - ■ -+X% +1 ) — 1, R(e) fc+1 ) above 
the points represented by .Mo using [2J Algorithm 12.18 (Parametrized 
Bounded Algebraic Sampling)]. 
Step 2. Call Algorithm [6] (Baby-giant Roadmap for Bounded Algebraic Sets) 
with input Q e , Mo, m = 0, performing arithmetic operations in the 
domain D[e]. The algorithm outputs a roadmap 

BGRM(Zer(Q e , R(e) k+l ),M ) 

composed of points and curves whose description involves e. 
Step 3. Denote by L the set of polynomials in D[e] whose signs have been 
determined in the preceding computation and take 

a = min c(P) 

(Notation 17. 2p . Replace e by a in the polynomial Q e to get a polyno- 
mial Q a . Replace e by a in the output roadmap to obtain a roadmap 
BGRM(Zer(Q a ,R fc+1 ),M J/i0 ). When projected to R k , this roadmap 
gives a roadmap 

BGRM(Zer(Q, R k ),M yfi ) n B k (0, 1 /a). 

Step 4. In order to extend the roadmap outside the ball B(0, 1/a) collect all 
the points (yi, . . . ,yk, Vk+i) £ R(e) fc+1 in the roadmap 

BGRM(Zer(Q e , R(e) k+1 ),M ) 

which satisfies e(y\ + . . . + yf.) = 1. Each such point is described by a 
real univariate representation involving e. Add to the roadmap the 
curve segment obtained by first forgetting the last coordinate and 
then treating e as a parameter which varies vary over (0, a] to get a 
roadmap BGRM(Zer(Q, R fc ), Mo). 

PROOF of correctness. The proof of correctness follows from the proof 
of correctness of Algorithm [6] (Baby-giant Roadmap for Bounded Algebraic 
Sets). □ 



Complexity analysis. The complexity is dominated by the complexity of 
Step 2. □ 
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Proofs of Theorem and Corollary \1.S\ . Follows directly from the correct- 
ness and complexity analysis of Algorithm [7] (Baby-giant Roadmap for Gen- 
eral Algebraic Sets) , after substituting m = and p = \[k. □ 

8. Appendix: computing the limit of bounded points and curve 

segments 

8.1. Limit of a bounded point. Before computing the limit of a bounded 
point we need to explain how to perform some useful computations modulo 
a quasi- monic triangular Thorn encoding J 7 , a representing a point t £ R m . 
We associate to t € R m specified by a triangular Thorn encoding J 7 , a, 

F = (/[l]> ■ • • >0[m])> f[i] 6 D[Ti, . . . ,Tj], 

the ordered domain D[i] contained in R and generated by t. 

We now aim at describing the pseudo-inversion of a non-zero element in 
the domain D[t] specified by J 7 , a. 

Definition 8.1. A pseudo-inverse of / G T)[t] is an element g € D[i] such 
that fg £ D is strictly positive. 

This notion is delicate as the computation of the pseudo-inverse sometimes 
requires us to update the triangular Thorn encoding specifying t, in the 
spirit of dynamical methods in algebra (see for example [8]). We start with 
a motivating example. 

Example 8.2. We consider t, specified as the root of 

f(T) = T 4 - T 2 - 2 

giving signs (+, +, +, +) to the set Der(/) of derivatives of /. 

Consider T 2 + 1. It is easy to see, using for example [21 Algorithm 10.13 
(Sign Determination Algorithm)] applied to / and the list Der(/),T 2 + 1, 
that the sign of T 2 + 1 at t is positive. In order to compute its pseudo- 
inverse, we perform [21 Algorithm 8.22 (Extended Signed Subresultant)] of 
/ and T + 1. If f(T) and T 2 + 1 were coprime, we would obtain the 
pseudo-inverse of T 2 + 1 modulo f(T) since the last subresultant would be 
a non zero constant in D. But f(T) and T 2 + 1 are not coprime and their 
gcd is T 2 + 1. bo we divide f(T) by T 2 + 1, obtain a new polynomial 
g(T) = T 2 — 2 and check that the root t of f(T) giving signs (+,+,+,+) 
to the set Der(/) coincides with \[2 which is the root of T 2 — 2 making 
the derivative g'(T) = 2T positive, using again -for example- [21 Algorithm 
10.13 (Sign Determination Algorithm)]. It is now possible to pseudo-reduce 
T 2 + 1 modulo g(T), which gives 3. 

In other words, during the process of computing the pseudo-inverse of 
T 2 + 1 we discovered the factor g(T) of f(T) having t as a root and coprime 
with T 2 + 1. Using this new description of t we have been able to compute 
a pseudo- inverse of T 2 + 1. 

We can now describe the computation of the pseudo-inverse in general. 
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Description 8.3. Given t = (t\,..., t m ) £ R m specified by the quasi-monic 
triangular Thorn encoding T = (/[i], • • • , f[ m ]), cr = (f 1, • • • , Cm), we describe 
how to compute a pseudo- inverse of a non-zero elements of D [t] . 

We proceed by induction on the number m of variables of J- . 

If m = there is nothing to do since D is an ordered domain. 

If m ^ 0, let t' = (ti, . . . ,t m -i) specified by T = {f\\\,. . . ,/[ m _i]),cr = 
(<7l, . . . ,cr m _i). 

We consider / as a polynomial in T m whose coefficients, which are ele- 
ments of 

{h G D[Ti, . . . , T m _i] | deg Ti (fc) <deg T .(/ [i] ),i = l,...,m-l} 

represent elements of D[i']. 

We first decide the sign of / at t, which is done by [21 Algorithm 12.19 
(Triangular Sign Determination Algorithm)]. 

If f{t) 7^ 0, we try to pseudo-invert / modulo T. We perform [21 
Algorithm 8.22 (Extended Signed Subresultant)] for / and /[ m ], with re- 
spect to the variable T m and compute a gcd(/, /r m i) € D[i'] (the last non 
zero subresultant polynomial) as well as the cofactors u, v £ D[t'] with 
uf + vf [m] = gcd(/,/ [m] ). 

(1) If gcd(/, /r m i) is of degree in T m , u is a quasi-inverse of /. 

(2) If gcd(/, /r m i) is of degree > in T m , we have discovered a factor 
of f\ m ] ■ We define h as the quasi-monic polynomial proportional to 
/[■m]/ gcd(/, /[ m ]) obtained by Algorithm 8.22 (Extended Signed 
Subresultant)] (see [H Algorithm 10.1 (Gcd and Gcd-free part)]). 
We perform [21 Algorithm 12.19 (Triangular Sign Determination)] 
applied to /r m ] and Der(/r m i), Der(/t) to identify the Thorn encoding 
r of ir m ] as a root of h. We replace /r m i by /i and <7[ m ] by r in J 7 . Now 
/ and the new /r m i, considered as polynomials in Tr m ] are coprime 
and we can invert / module /r m i . 

Proposition 8.4. Let D 6e an ordered domain contained in a real closed 
field R, and t = (t%, . . . ,t m ) S R m specified by a triangular Thorn encoding 
J 7 , a, 

F = (/[i] i • • • i f[m] ) > f[i] G D[Ti , . . . , Tj] . 
Let d be a bound of the degree of f^ with respect to each Tj, i = 1, . . . , m, j < 

i. 

a) Lf f £ D[Ti, . . . ,T m ] is a polynomial of degree D, the complexity of 
computing the pseudo-reduction PsRed(/, F), is (0((D — d)d)) m arithmetic 
operations in D. 

b) The complexity of the computation of the pseudo-inverse of an element 
of~D[t] is d 0<yTri ^ arithmetic operations in D. 

Proof. The proof of a) is made by induction on the number of variables m 
of T . If m = 0, there is nothing to do. If m / 0, the pseudo-division 
process by g^ m j takes 0((D — d)d) arithmetic operations in D[Ti,...,T m ] 
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using the complexity analysis of [21 Algorithm 8.3 (Euclidean Division)] and 
produces d polynomials of degree D + (D — d + l)d in D[Ti, . . . , T m ], so the 
complexity is {0{{D — d)d)) m . The end of the reduction of / coincides with 
the reduction of these d polynomials modulo 

? = (/[1]>- • • ) /[m-l])j f[i] G D [Ti, • • • > T i], 

and costs d(0((D — d)d))^ m ~ 1 ^ by the induction hypothesis. 

The proof of b) proceeds also by induction on the number of variables m 
of T. 

If m = 1, the computation of a gcd takes (d+ l) c operations in the domain 
D, for some universal constant c > 0, using the complexity analysis of [21 
Algorithm 8.22 (Extended Signed Subresultant)] and [21 Algorithm 10.13 
(Sign Determination)]. 

If m > 1, let i = (t', «), and we suppose by induction hypothesis that the 
complexity of arithmetic operations in D[t'] is (d + l) c ( m_1 ) arithmetic oper- 
ations in the ordered domain D. The statement is clear since the arithmetic 
operations in the domain ~D[t] are using (d + l) c operations in the domain 
D[i'] using the complexity analysis of [21 Algorithm 8.22 (Extended Signed 
Subresultant)] and [21 Algorithm 10.13 (Sign Determination)]. □ 

We can now give the description of Algorithm [3] (Limit of a Bounded 
Point). 

Procedure. Remove from g e (T,U) the monomials vanishing at t, using [21 Al- 
gorithm 12.19 (Triangular Sign Determination)]. Supposing with- 
out loss of generality that all the coefficients of g e (t,U) are not 
multiple of e, denote by g{T\, . . . ,T m , U) the polynomial obtained 
by substituting to e in g e {T\, . . . ,T m ,U). Similarly denote by 
G(T±, . . . ,T m , U) the polynomials obtained by substituting to e 
in G s (Ti, . . . ,T m , U). 

Compute the set S of Thorn encodings of roots of g(t, U) using [21 
Algorithm 12.19 (Triangular Sign Determination)]. Denoting by \i a 
the multiplicity of the root of g(t, U) with Thorn encoding a, define 
G a as the \i a — 1-th derivative of G with respect to U. 

Identify the Thorn encoding a and G a representing z using [21 Al- 
gorithm 12.19 (Triangular Sign Determination)], by checking whether 
a ball of infinitesimal radius 5 (1 S> S S> e > 0) around the point 
x represented by the real univariate representation g, a, G a contains 

Pseudo-invert the leading coefficient of the univariate representa- 
tion, denote by J 7 ', a 1 the new triangular Thorn encoding describing 
t and pseudo-reduce by T' . 



Complexity analysis: Follows from the complexity of [21 Algorithm 12.19 
(Triangular Sign Determination)]. □ 
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8.2. Limit of a curve. Computing the limit of a curve is not immediate 
when some part of the curve has a vertical limit, as seen in the following 
example. 

Example 8.5. Consider the semi-algebraic curve 7 : [0, e] —> R(e) 3 , parametrized 
by the X\ coordinate defined by 

7(37) = (xi, 72 Oi), 73 (xi)), x\ e [0,e] 

where (72(21), 73(27))) is the solution of the triangular system, 

X 2 - xx/e = 0, 
X\ + Xi - 1 = 0, 

with Thorn encoding (0, +), (0, +, +). 

Notice that the image of 7 is contained in the cylinder of unit radius with 
axis the Xi-axis and is bounded over R. The image of 7 under the lim £ map 
is contained in a circle in the plane Xx = 0, and can no longer be described 
as a curve parametrized by the Xi-coordinate. 

However, it is possible to reparametrize 7 by the ^-coordinate. By 
doing so we obtain another semi-algebraic curve <p : [0, 1] -»• R(e) 3 (having 
the same image as 7) defined by 

lf(x 2 ) = (ipx(x2),X2,¥3(x2)),X 2 G [0,1] 

where {ipx{x 2 ), 933(22)) is the real solution of the triangular system 

Xx - ex 2 = 0, 
X$ + x%-l = 0, 

with Thorn encoding (0, — ), (0, +, +). Notice that the image under lim e 
of the curve which is the graph of <p can be easily described as the curve 
represented by the following triangular system parametrized by x 2 G [0, 1] 

Xx = Q, 
Xi + x 2 2 -l = 0, 

and Thorn encoding (0, —1), (0, +, +). 

This is the reason why some kind of reparametrization is necessary before 
computing the limit. 

8.2.1. Reparametrization of curve segments. We define the notion of well- 
parametrized curve, and prove that the limit of a well-parametrized curve 
is easy to describe. 

Definition 8.6. A differentiable semi-algebraic curve 



7 = ( 7 i,...,7fc) : (a,b)^R k 
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parametrized by X\ (i.e. 71 (x±) = x±) is well-parametrized if for every 
x\ G (a, 6), 

SO*" 

Let t G R m be represented by a triangular Thorn encoding T , <r, and 

be a curve segment with parameter Xj over t on (01,02) where 01 and 02 
are the elements of R represented by the Thorn encodings fx, o\ and /2, o~2. 
The curve segment 

is well-parametrized if the semi-algebraic curve 7 : (01,02) — >• R(e) fc defined 
by 

/ x = / 9i(t,Xj,u{xj)) g k (t,Xj,u(xj)) 
3 \9o(t, Xj,u(xj)) ' " ' g (*, Xj,u(xj)) 

is well-parametrized, where it : (01,02) — > R fc_( ^+ r ) maps each Xj G (01,02) 
to the root of git, xj, U) with Thorn encoding r. This means that 

-yff gi(*,sjXsj)) y\ 2 < fe 

where the derivative is taken with respect to xj. 

Example 18.51 is not a well-parametrized curve segment. 

If a curve segment defined over R(e) is well-parametrized, and represents 
a curve bounded over R, then the image of the curve under the lim £ map 
can be easily described. The following proposition explains why this is true. 

Proposition 8.7. Let (a £ ,b £ ) C R(e), r < j < k, z £ = (z £t x, ■ ■ ■ , z Er ) G 
R(e) r , and 

le = (Ze,le,r+1, ■ ■ ■ ,le,k) ■ (o e , 6 e ) -> {z e } X R(e) k ~ r 

a semi- algebraic differentiate curve parametrized by Xj and bounded over 
R. If 7 E is well-parametrized, 

(1) /or eac/t x G (lim £ a £ , lim £ b £ ) and any x e G (a £ ,b £ ) with lim e x £ = x, 
7(3;) := lim £ 7 £ (x) = lim £ 7 £ (x £ ), 

(2) lim £ 7 e ((a £ ,6 £ )) = 7([lim £ a £ , lim £ b e ]). 

In other words, the graph of the semi-algebraic function 7(3:) := lim £ 7 e (x) 
is the image under lim £ of the graph of 7 £ . 

Proof. It follows from the definition of being well-parametrized that \ \j' £ (x) \ \ < 
y/k for all x G (a £ ,b £ ). By the semi-algebraic mean value theorem [2j Ex- 
ercice 3.4] we have that for each x G (lim e a £ , lim e b £ ) and any x £ G (a £ ,b £ ) 
with lim £ x £ = x, 

\\le(x) -7 e (x e )|| = ||7 £ (w £ )|||x-x £ |, 
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for some w G (x,x e ) (assuming without loss of generality that x < x e ). 
Taking the image under lim e and noticing that ||7e(w e )|| is bounded over 
R by the previous observation, we obtain that lim £ 7 £ (x) = lim e 7 £ (x £ ), 
proving (1). This implies that the function 7 : (lim £ a £ , lim £ b £ ) — > R fc defined 
by j(x) = lim e 7 e (x) is a continuous, bounded (since 7 £ is bounded over 
R) semi-algebraic function, and hence can be extended to a continuous, 
bounded semi-algebraic function on the closed interval [lim e a E , lim £ b e ], and 
(2) follows. □ 

A semi-algebraic curve is in general not well-parametrized. However, 
subdividing if necessary the curve into several pieces, it is possible to choose 
for each such piece a parametrizing coordinate which makes the piece well- 
parametrized. This is what we do in Algorithm (Reparametrization of a 
Curve) . 

Algorithm 8. [Reparametrization of a Curve] 

Input. (1) t £ R m represented by a triangular Thorn encoding J-, a, 
(2) a bounded curve S represented by a curve segment, 

fl,<ri,f2,<T2,9,T,G 

with parameter X\ in R fc over t on (a, b). 
All the polynomials in the input have coefficients in D. 
Output. (1) A finite set V = \y\, . . . , vn-i}, of real univariate representation 
over (t, q), where each a is a Thorn encoding over t fixing X m a) . 
(2) A finite set W = {wi, . . . , iojv}, of curve segments with wi 
parametrized by . 

Moreover, the union of the curves represented by W, and the 
points represented by V define a partition of S. 
Complexity If the polynomials occurring in the input have degrees bounded by 
D then the complexity of the algorithm is bounded by D°^ m \ 

Procedure. 

Step 1. Let gi(Xi,T) = Xigo(Xx,T), and for each i, 1 < i < k, let 
( dg\ ( dgi dg \ ( dg { dg \ ( dg 



dT J \dXi yt dXj V 5T Ul dT J \dX x 

(which is proportional to the projection on the i-th coordinate of the 
tangent vector to the input curve by the chain rule) and 

i=i 

Step 2. Computing REliniT(G ; i, 5), 1 < i < k, using [21 Algorithm 11.19 (Re- 
stricted Elimination)], obtain a family C of polynomials in D[Ti, . . . , T„ 
Subdivide (a, b) in a finite union of points and intervals over which 
the signs of the polynomials in C are fixed using [21 Algorithm 12.23 
(Triangular Sampling Points)] and get a = c\ < . . . < cl = b, where 
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each Cj is represented by a Thorn encoding (Cj,aj) over t £ R m , 
such that for each j, 1 < j < L, there exists an £(j), 1 < t(j) < k, 
such that for all x\ E (cj_i,Cj), G^{t, x\, u{x{)) > 0, denoting by 
u(x\) the root of g(t,xi,U) with Thorn encoding r. For each j fix 
an £(j) satisfying this property. 
Step 3. For each j, 1 < j < L, reparametrize the segment of the input curve 
over the interval (cj_i, Cj) using the coordinate X^y Suppose with- 
out loss of generality from here on that £(j) = 2. 
Step 3 a). Set 

H : = g 2 + (X 2 ■ g (T, X l , U) - g 2 (T, X\, U)) 2 € D [T, X U X 2 ,U]. 

Note that Zer(H(t, -),R 3 ) is a curve bounded over R (by assumption 
on the input). Call Algorithm [1] (Curve Segments) with input the 
polynomial H, and the triangular system J 7 , a, noticing that X 2 is 
now the parameter. 
Step 3 b). For each element 

(h{T,X 2 ,V),a h ,H(T,X 2 ,V)) G V u 

where 

H(T, X 2 ,V) = (h (T, X 2 ,V), hi (T, X 2 , V), h 2 (T, X 2 , V)), 

and T>i is a set of distinguished points in the output of the call to 
Algorithm [1] (Curve Segments) in the previous step, use [21 Algo- 
rithm 12.19 (Triangular Sign Determination)] to check if the point 
(xi,x 2 ,u) represented by (h,ah,(ho,hi,X 2 ho,h 2 )) over i, coincides 
with 

goi^x^uix!)) 
Retain only the element 

(h(T,X 2 ,V),a h ,H(T,X 2 ,V)) €Vt 

for which this is the case, and add to the set V the real univari- 
ate representation u = (h, er^, Gh) (see Notation I4.6|) representing a 
point Vh G R fc , with parameter X 2 over t. 
Step 3 c). For each element 

Ux{T,V),a x J 2 (T,V),G 2 ,h{T,X 2 ,V),a h ,H(T,X 2 ,V))eC h 

where 

H{T,X 2 , V) = (h (T, X 2 ,V), hi(T, X 2 ,V), h 2 (T, X 2 ,V)) 

use [21 Algorithm 12.19 (Triangular Sign Determination)] to check if 
the point (x\,x 2 ,u) represented by 

h(T,X 2 ,V),cr h , (h ,hi,X 2 h ,h 2 ) 
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over t, coincides with 

g 2 {t,x 1 ,u{xi)) 

[XI, —77 7 vT,M(xi)) 

for X2 = (vi + v 2 )/2 where v\,v 2 are represented by (/i,ci) and 
(/2)C2) respectively. Retain only the element of Ci for which this is 
the case, and add to the set W the curve segment, or, f 2 ,a 2 , h, ah, Gh), 
with parameter X2 over t (see Notation 14. 6p . 

PROOF of correctness. Let {fi,<Ji, f 2 ,a 2 ,g,r,G) be a curve segment 
parametrized by X\ over t representing the curve 7 : (a, 6) — >• R fc . 
Let (c, d) be a sub-interval of (a, 6) such that for every x\ G (a, 6) 

(8.1) Gt(t, x x , u{x 1 )) = kFf{t, xi, u(xi)) - if (i, xi, «(xi)) > 0. 

j=i 

(using the notation of Step 1 and Step 2). 
Since 



axi 



1 

~7F 



the mapping 7^ from (c, d) to (c', </); with c' = Ji(c), d! = ji(d) is invertible. 
Defining j(xe) = 7(7^ (%t)), 7((c',d')) = 7(( c i^)) i s well-parametrized by 

Moreover, at each point x\ G (a, b) such a choice of £ exists, since there 

( \ 

must exist an £, 1 < £ < k such that ( ^ v ) is at least the average value 



In Step 2 of the algorithm we obtain a partition of the 
k t— f V / 

interval (a, 6) into points and open intervals, such that over each sub-interval 
(cj-i,Cj) of the partition, there exists an index £ = £(j) such that (|8.1j) is 
satisfied at each pointu G (cj_i,Cj). 

Each curve segment corresponding to elements of V output by the algo- 
rithm is thus well-parametrized. The remaining property of the output is 
a consequence of the correctness of Algorithm Q] (Curve Segments), and [21 
Algorithm 12.19 (Triangular Sign Determination)]. □ 

Complexity analysis. 

Let D be a bound on the degrees of the polynomials in the input. The 
complexity of Steps 1 and 2 is bounded by k°^ D°( m ) from the complexity 
of [21 Algorithm 11.19 (Restricted Elimination)], and [21 Algorithm 12.23 
(Triangular Sampling Points)], noting that the number of polynomials in C 
is bounded by k°^D°( m \ 

In Steps 3-4 the Algorithm [T] (Curve Segments) and [21 Algorithm 12.19 
(Triangular Sign Determination)] are both called with a constant number 
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of variables in the input. Using the complexity analysis of these algorithms, 
the total complexity is bounded by D°( m ' . □ 

8.2.2. Limit of a curve. We are now ready to describe Algorithm U] (Limit 
of a Curve). The algorithm proceeds by reparametrizing the curve and 
computing the limit of the well-parametrized curve segments so obtained, 
as explained below. 

Procedure. 

Step 1. Let T = (T±, . . . ,T m ), X' = (X\, . . . ,X r ). Call a slight variant of 
[2j Algorithm 12.18 (Parametrized Bounded Algebraic Sampling)], 
pseudo-reducing intermediate computations modulo T using (|4.ip . 
with input 

£ A 2 € D[e, T, X'] 

and parameters e, T, and output the set IA £ of parametrized univari- 
ate representations with variable U. 

For every (h £ ,H £ ) £ U £ , use [2j Algorithm 12.20 (Triangular Thorn 
Encoding)] with input the triangular system (J 7 , h £ ) to compute the 
Thorn encodings of the real roots of h e (y, U). If 

%e = ■ • • i h[ r ]) 

with hfi] 6 D[T,X\, . . . ,Xi,] substitute the variables X' in 

U Der ^(%i) 

l,...,r 

using H £ by fj4.5j) and define a family A of polynomials in e, T, ?7. 
Using j2j Algorithm 12. (Triangular Sign Determination)], compute 
the signs of the polynomials of A at the roots of h £ (y, U). Comparing 
the Thorn encodings, identify a specific (h £ ,r s , H e ) representing z £ 
over t. 

Then apply Algorithm [3] (Limit of a Bounded Point) with input 
(he,T £ ,H £ ) representing z £ over t to obtain a real univariate repre- 
sentation p z , p z , P z representing z over t. 

Step 2. Using Algorithm [8] (Reparametrization of a Curve) reparametrize 
the input curve segment. 

Step 3. For every well parametrized curve segment output in Step 2, S E , 
represented by 

(/e,l) °~e,l, fe,2, C e ,2, 9e, T £ , G £ ), 

do the following. 

First reorder the variables to ensure that the parameter of S £ is 
X r +i- 

Then compute a description of lim £ (S' £ ). This process is going to 
generate a finite list of open intervals and points above which the 
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representation of the restriction of the curve lim £ (5 e ) by a curve 
segment is fixed. This is done as follows. 
Step 3 a). Denote by a £t \ the element of R(e) represented by / e ,i(e, T, X' , X r+ i),a £t \ 
over (t,z e ). 

Call a slight variant of Algorithm 12.18 (Parametrized Bounded 
Algebraic Sampling)], pseudo-reducing intermediate computations 
modulo J- using (14. ip . with input 

A 2 + f £ ,i{T, X',X r+1 ) 2 G D[e, T, X' , X r+1 ] 

and parameters e, T, and output a set £/ £ of parametrized univariate 
representations with variable U. 

For every (h £ ,H £ ) G Z^, use [21 Algorithm 12.20 (Triangular Thorn 
Encoding)] with input the triangular system (J 7 , h £ ) to compute the 
Thorn encodings of the real roots of h £ (y, U). 

If 

7-L £ = (h[i],. . . , h[ r ]) 
with hu] G D[T, Xi, . . . , Xj, ] substitute the variables X' , X r+ i in 

Der Xr+1 (/ £il ( £ ,T,^X r+1 ),X r+1 U J Der X( (fc w ) 

l,...,r 

using (|4.5p and define a family i3 of polynomials in e, T, [/. Using [21 
Algorithm 12. (Triangular Sign Determination)], compute the signs 
of the polynomials of B at the roots of h e (y,U). Comparing the 
Thorn encodings, identify a specific (h £ ,T £ , H £ ) representing (z £ , a e ,i) 
over t. 

Then apply Algorithm [3] (Limit of a Bounded Point) with input 
(h e , t £ ,H £ ) representing (z e , a e ,i) over t to obtain a quasi- monic real 
univariate representation p z , ai , Pz,a x , Pz,a^ representing (z,a±) over 
t with qi = lim e (a £| i). Obtain a Thorn encoding over t, of «i using 
[21 Algorithm 15.1 (Projection)]. 

Similarly, for a £ ^ the element of R(e) represented by f e ,2(T, X' , X r+ i), cr £i 2 
over (t, z e ), compute a Thorn encoding over t, of «2 = h m e(Q!e,2)- 
Step 3 b). Perform a slight variant of [21 Algorithm 12.18 (Parametrized Bounded 
Algebraic Sampling)], pseudo-reducing intermediate computations 
modulo J- using (I4.1D . with input 

Y, A * + 9e{T, X', X r+U V) 2 € D [e , T, X', A r+1 , E7] 

with parameters e, T, X r+ \ and output a set V £ of parametrized uni- 
variate representations with parameter e, T, X r+ \ and variable V. 
Denote by T £ the set of polynomials f £ such that there exists F £ 
such that (/ e ,F B ) G V E . Note that / E G D[e, T, X r+ i, V]. 
Step 3 c). Compute the family of coefficients C C D[T, A r+ i] of the polynomi- 
als f £ G T £ considered as elements of D[T, X r+ i][e, V] and the list 
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C of non-empty conditions = 0, 7^ satisfied by C in R using [21 
Algorithm 12.23 (Triangular Sample Points)]. Note that for every 
x r +i in the realization of r E C, the orders in e of the coefficients of 
the polynomials in T £ {t,x r+ i) C D[V] are fixed. For every f e £ T £ 
we denote by o(/ e ,r) the maximal order in e of the coefficients of 
f £ (t, av+i) on the realization of r and by T T C D[T, X r+ i, V] the set 
of polynomials obtained by substituting for e in £ - °u e > T ) f e . 
Step 3 d). Define 

T= \jT T cD[T,X r+1 ,V}. 

Compute 

£ = C\J RElirm/(/,Der(/)) c D[T,X r+1 ]. 

using [21 Algorithm 11.19 (Restricted Elimination)], so that the 
Thorn encodings of the real roots of f{t,x r+ \,V) are fixed when 
x r +i varies in an open interval defined by the roots of the polyno- 
mials £ (i). 

Step 3 e). Compute using [2j Algorithm 12.19 (Triangular Sign Determination)] 
the Thorn encodings of the real roots of the polynomials in £(t), and 
the ordered list c\ < ■ ■ ■ < Ch-i of the roots of the polynomials in 
£(t) in the interval (co,Ch), with Co = = «2- Denote by Cj,pj 

a polynomial in £ (t) and a Thorn encoding representing Cj. 

Step 3 f). For every j from 1 to h — 1 and for every / € J 7 , determine, using 
O Algorithm 12.19 (Triangular Sign Determination)] the Thorn en- 
coding f(t,Cj,V),Tj of a root Vj such that i>j = lim E (t> E ), where v £ 
is the root of / e (e, i, Cj, V) with Thorn encoding r e . The multiplicity 
/ij of the root Vj is determined by tj . 

Step 3 g). For every j from 1 to h, define I = (c^-i, c— j). For every / € J 7 de- 
termine, using [21 Algorithm 12.19 (Triangular Sign Determination)] 
the Thorn encoding //(i, x r+ i, V), r/ of a root vj(x r+ i), of multiplic- 
ity ^7 such that for every x r+ \ G /, t>/(x r+ i) = lim e (« e ) where u e is 
the root of f e (e, t, x r+ i, V) with Thorn encoding t £ . The multiplicity 
Hi of the root vj(x r+ i) is determined by r/. 

Step 3 h) . Given (f e , F £ ) in U £ denote by (gp e ,Gf s ) the k — r + 1-tuple of poly- 
nomials obtained by substituting in (g £ ,G £ ) the variables X',U by 
F £ (see Notation [33]). Denote by V £ C D[e,T,Xj, V] the set of 
A; — r + 1-tuples of polynomials (gF s , Gf e )- 

Step 3 i). For every j from 1 to h — 1 and every (h £ ,H £ ) S V e , with ff e = 
(/i £j o, /i £ ,r+2 5 • • • j h £: k) determine the order in e of 

h E (e,t,Cj,Vj),h £i i(e,t,Cj,Vj). 

This is done by determining the signs of the coefficient he, hi^e of £ in 
h(e, t, Cj,vj), hi(e, t, Cj,Vj) using [21 Algorithm 12.19 (Triangular Sign 
Determination)]. Retain those (h £ ,H £ ) such that o(/i £i o) > o(h e ^) 
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for all i from r + 2 to k and replace e by in (e~ o( - h ^h £ , e~°^ h£fi ^ H £ ), 
which defines a set Tij. Inspecting every (h,H) G Hj, determine, 
using [21 Algorithm 12.19 (Triangular Sign Determination)], a k — 
r + 1-tuple (hj, Hj) with the following property. Let dj be the point 
represented by the real univariate representation 

{hj (T, X r+1 , V) , Tj , H^- X) (T, X r+1 , V)) 

over t, u. The image under lim e of the point of S £ with X T+ \- 
coordinate (cj) is (z,Cj,dj). 
Step 3 j). For every j from 1 to h define I = (cj-i, c — j). For every (h £ , H £ ) G 
Vg, with iJ £ = (/i E| o, h £)r+ 2, . . . , /i 6 ft) subdivide 7 so that the order in 
£ of h £ (e,t,Cj,Vj) and /&j(e,i, x r +i,«j(ai r +i)) is fixed. This is done 
by computing 

£i= (J RElim v (/, c D[T, X r+1 ], 

/eJ r ,(ft,H e )eV^,o<^<deg e ^ 

and 

f j,i = [J RElirm/(/, ^) C D[T, X r+1 ], 

feJ 7 ,{h,H s )eV e ,0<£<deg e hi 

using [21 Algorithm 11.19 (Restricted Elimination)]. 
Defining 

£'i = £1 [J 

ie0,r+2,...,k 

compute the Thorn encodings of the roots of the polynomials in 
£j{t), using [21 Algorithm 12.19 (Triangular Sign Determination)]. 
On each open interval J between two successive roots, the order in 
e, denoted by o(h £ ),o(h £ ^) of the polynomials 

h(e, t, x r+1 ,vj(x r+ i)), hi(e, t, x r+1 ,vj(x r+1 )) 

remains fixed. Retain those (h £ ,H £ ) such that o(h £ $) > o{h £ ^) for 
all i from r + 2 to k and replace e by in e~°( hefi \h, H £ ), which 
defines a set %j. Inspecting every (h,H) G %j, determine, using (2[ 
Algorithm 12.19 (Triangular Sign Determination)], a k — r + 1-tuple 
(hj,Hj) such that the point represented by 

{hj(t, x r+ i,v I ),H < f J ~ 1 \t, Xr+i, u/)) 

is the image under lim £ of the point of S £ with X r+ i-coordinate x r +i, 
where fij is the multiplicity of uj{x r+ \) as a root of hj(x r+ \, V). 

Let wj be the curve represented by the curve segment represen- 
tation 

hi (T, X r+1 , U) , Tj , ^ J_1) (T, X r+1 , CT) 
with parameter X r+ \ over i, tt. 
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Step 3 k). Let c\ < ■ ■ ■ < cn-i denote the set of all the elements of R computed 
in Steps 2 d), and 2 i) above, and cn = c. Re-index each Vj computed 
in Step 3 h), such that dj lies above Cj. Similarly, re- index each wi 
computed in Step 3 i) by some j, 1 < j < N, so that Wj lies above 
the interval (cj_i,Cj). 

Output the list of di, . . . , d^-i, and w±, . . . , wn- 

Proof of correctness. Let <y £ : (a £i i,a £i2 ) — > R(e) k be the curve repre- 
sented by a well-parametrized curve segment 

/e,l> ^T' / £j2 , C £ ,2, g e ,T £ , G £ 

computed in Step 2. 

Let G : (a±, «2) — s- R fc be the curve whose image equals the image of r y £ 
under lim e . Since the input curve segment is well-parametrized it follows 
from Proposition 18.71 that in order to compute for any x\ £ (co, cat), G{x\) it 
suffices to compute lim £ 7 e (xi). The proof of correctness of the algorithm is 
then similar to the proof of correctness of Algorithm [3] (Limit of a Bounded 
Point). 

□ 

Complexity analysis. Let D be a bound on the degrees of all polynomials 
appearing in the input. We first bound the degrees in the various variables, 
e, T, X' , X r+ \, U, V of the polynomials computed in various steps of the al- 
gorithm. In Step 1, the degrees of the polynomials in U £ are bounded as 
follows. The degrees in e, U are bounded by by the complexity analy- 

sis of [21 Algorithm 12.18 (Parametrized Bounded Algebraic Sampling)] and 
the degrees in the Tj are bounded by D, because of the pseudo-reduction. 
Moreover, the complexity of this step is bounded by £)°( m + r ) from the com- 
plexity of [21 Algorithm 12.18 (Parametrized Bounded Algebraic Sampling)] 
and the complexity of pseudo-reduction (see Definition 14. 2|) . 

The degrees in e,Ti,X',U in the output of Step 2 are all bounded by 
D°W and the complexity of Step 2 is bounded by 

ih - r )°( 1 ) £)0(m+r) = k O(l) D 0{m+r) 

using the complexity analysis of Algorithm [8] (Reparametrization of a Curve). 

The degrees of the polynomials in Step 3 a are bounded as follows. In the 
output of the call to [21 Algorithm 12.18 (Parametrized Bounded Algebraic 
Sampling)], the degrees in e,U are bounded by D 0< - r \ and the degrees in 
the Tj are bounded by D. Now, from the complexity analysis of Algorithm 
[3] (Limit of a Bounded Point) it follows that the degrees in the Tj of the 
polynomials output are bounded by D and those in e, U are bounded by 
D°( r \ Moreover, the complexity of Step 3 a is bounded by Jj°( m + r ) from 
the complexity of [21 Algorithm 12.18 (Parametrized Bounded Algebraic 
Sampling)], the complexity of Algorithm [3] (Limit of a Bounded Point) and 
the complexity of pseudo-reduction (see Proposition I8.4h . 
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The degrees of the polynomials in Step 3 b are bounded as follows. In the 
output of the call to [2j Algorithm 12.18 (Parametrized Bounded Algebraic 
Sampling)], the degrees in e,X r+ \,V are bounded by D°( r \ and the de- 
grees in the Tj are bounded by D. The complexity of Step 3 b is bounded by 
£)0(m+r) f rom complexity of [2j Algorithm 12.18 (Parametrized Bounded 
Algebraic Sampling)], and the complexity of pseudo-reduction (see Defini- 
tion 

The complexity of Step 3 c is bounded by r)°( m + r ) using the degree 
bounds from the complexity analysis of the previous steps and the complex- 
ity of [21 Algorithm 12.23 (Triangular Sample Points)]. 

It now follows from the complexity analysis of [21 Algorithm 12.19 (Trian- 
gular Sign Determination)], [21 Algorithm 11.19 (Restricted Elimination)], 
and the degree estimates proved above that the complexity of the remain- 
ing steps are all bounded by k 0( ^D°^ m+r \ Thus, the complexity of the 
algorithm is bounded by D 0{ - m+r \ □ 
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